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Chapter  1 


INTRODUCTION 


Atherosclerosis  is  a  condition  of  surface  degradation  of  the  lining  of  the  ar¬ 
teries.  This  condition  can  lead  to  heart  attack  or  stroke.  Heart  attacks  alone 
cause  500,000  deaths  a  year  in  the  United  States.  The  detection  of  atheroscle¬ 
rosis  early  in  its  development  can  improve  the  prognosis  of  the  individual  with 
the  condition.  The  use  of  ultrasound  is  one  method  to  detect  atherosclerosis. 
The  thesis  consists  of  six  parts.  The  first  part  describes  the  disease  and  cur¬ 
rent  ultrasound  modalities  used  for  the  diagnosis  of  the  disease.  The  second 
part  discusses  the  theoretical  basis  of  scattering  of  waves  from  rough  surfaces 
and  how  it  can  be  applied  to  studying  atherosclerosis.  The  third  part  discusses 
the  characteristics  of  phantoms  and  tissues  that  were  experimentally  studied. 
The  fourth  section  discusses  the  equipment  and  procedures  used.  The  fifth 
section  presents  the  results.  The  sixth  section  contains  the  conclusions.  The 
appendices  contain  supporting  data  and  the  software  used  for  collecting  and 
analyzing  data. 


Atherosclerosis 


To  discuss  the  disease  some  background  is  necessary  of  the  structure  of  an 
artery.  The  artery  is  composed  of  three  layers  of  tissue  (Gray,  1974,  p.  1126). 
The  innermost  layer,  the  tunica  intima,  itself  also  contains  three  parts.  The 


innermost  part  is  composed  of  smooth  epithelial  cells  oriented  longitudinally 
to  the  vessel.  These  epithelial  cells  serve  as  a  barrier  for  blood  components 
and  a  boundary  surface  that  has  little  drag  to  impede  the  blood  flow.  Also 
part  of  the  intima  is  a  subendothelial  layer  of  connective  tissue  with  branched 
cells  in  the  interspaces  of  the  tissue,  a  layer  which  is  more  developed  in  the 
larger  vessels.  The  chief  thickness  of  the  intima  is  due  to  a  third  part  which  is 
composed  of  elastic  longitudinal  fibers  with  small  elongated  perforations.  The 
next  layer  is  the  tunica  media.  The  media  is  composed  of  muscle  cells  that  are 
oriented  circumferentially  to  the  vessel.  These  muscles  control  the  diameter 
of  the  artery  depending  on  the  required  needs  of  the  organs  served  by  the 
artery.  The  outermost  layer  is  the  tunica  adventitia ,  composed  of  connective 
tissue  that  maintains  the  integrity  of  the  artery  and  its  maximum  diameter. 
This  layer  is  largely  absent  in  the  larger  arteries. 

The  exact  causes  of  the  disease  are  not  completely  defined,  although  some 
contributing  factors  are  known.  Chazov  and  Smirnov  documented  five  prob¬ 
lems  that  researchers  have  directed  their  efforts  at  solving  (1982,  p.  viii): 

1.  The  mechanism  regulating  cell  and  chemical  homeostasis. 

2.  Molecular  mechanisms  of  endothelial  injury  and  its  repair. 

3.  Mechanisms  of  intima  hyperplasia. 

4.  Lipidogenic  and  thrombogenic  mechanisms  of  plaque  formation. 

5.  Origin  of  foam  cells  and  pathways  of  extracellular  matrix  formation 
in  the  plaque. 

These  axe  all  basic  areas  for  research.  The  development  of  the  disease  is 
assumed  to  come  from  two  etiologies  (Likar  and  Robinson,  1985,  p.  12).  The 
first  is  insudation,  the  passage  and  accumulation  of  plasma  components  into 
the  intima.  The  second  etiology  is  encrustation,  the  accumulation  of  thrombi 
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or  blood  clots  at  points  of  arterial  weakness  or  damage  which  can  occur  from  a 
number  of  causes.  Atherosclerosis  may  result  in  plaques,  sometimes  occluding 
the  vessel,  containing  a  matrix  loaded  with  lipids  (Likar  and  Robinson,  1985, 
p.  10).  In  general  the  development  of  the  atherosclerotic  condition  begins 
with  a  thickening  of  the  intima.  This  is  accompanied  with  changes  in  the 
media,  either  local  or  diffuse  depending  on  external  factors  such  as  blood 
pressure,  thickening  rate,  surrounding  tissue  and  type  of  metabolic  disorder. 
This  increased  growth  is  not  accompanied  by  growth  of  the  vasa  vasorum, 
nutrient  vessels  supplying  the  larger  blood  vessels.  This  loss  of  nutrient  results 
in  regressive  metamorphosis  (Likar  and  Robinson,  1985,  p.  11).  Attempts 
to  repair  the  tissue  lead  to  clot  formation  and  further  thickening  with  the 
clots  having  a  fibrous  nature.  Ultimately,  the  vessel  becomes  partially  or 
totally  occluded,  also  accompanied  by  calcification  which  is  characterized  by 
granular  or  plate-like  shells  in  the  atheroma.  These  developmental  stages  are 
sometimes  used  as  references  for  characterization  of  diseased  tissue.  Barzilai  et 
al.  used  four  stages  of  development:  normal,  fibrous,  fibrofatty,  and  calcified 
(1987,  p.  460).  Picano  et  al.  used  five  stages:  normal,  fatty,  fibrofatty, 
fibrotic,  and  calcified  (1985,  p.  572). 

Another  more  recent  appraised  of  the  development  of  atherosclerosis  is 
found  in  a  1991  report  published  by  the  U.  S.  Department  of  Health  and  Hu¬ 
man  Services  (p.  9).  The  deposition  of  lipid  begins  in  the  arteries  of  children 
as  fatty  streaks.  The  extent  of  the  fatty  streaks  increases  in  adolescence  and 
early  adulthood.  Fibrous  plaques  (larger  lipid  deposits  with  a  covering  of 
smooth  muscle  and  connective  tissue)  adso  begin  to  appear  and  increase  dur¬ 
ing  the  second  and  third  decades.  Some  of  the  plaques  develop  other  changes 
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such  as  vascularization,  hemorrhage,  and  ulceration.  These  predispose  the 
vessels  to  thrombosis  during  middle  age  and  later. 


Imaging  of  Atherosclerosis 


Diagnosis  of  atherosclerosis  is  dependent  on  the  ability  to  image  the  lesions 
clinically.  There  are  both  noninvasive  and  invasive  intravascular  techniques  for 
imaging  atherosclerosis.  And  for  both  techniques  there  are  several  modalities, 
each  having  some  advantages  and  disadvantages. 

Noninvasive  techniques  include  contrast  angiography  and  B-mode  ultra¬ 
sound.  Angiography  involves  the  injection  of  a  radiopaque  dye  into  the  vas¬ 
cular  system  and  using  a  fluoroscope  or  x-ray  to  visualize  the  artery.  This  can 
be  real-time  with  fluoroscopy,  or  a  single  image  or  set  of  images  with  radiogra¬ 
phy.  The  technique  can  detect  stenosis,  or  narrowing  of  the  artery,  which  can 
be  measured  and  compared  with  adjacent  areas  of  the  vessel.  The  technique 
has  some  disadvantages:  it  can  not  distinguish  the  type  of  atherosclerotic  le¬ 
sion,  and  the  actual  measurement  may  not  be  the  maximum  stenosis  due  to 
orientation  of  the  image  plane,  as  most  lesions  are  eccentric  in  cross-section. 
This  can  however,  be  partially  corrected  by  orienting  the  patient  when  using 
fluoroscopy.  There  is  also  the  clinical  risk  associated  with  injection  of  contrast 
media  (Wells,  1982,  pp.  46-48). 

Imaging  of  atherosclerotic  plaque  is  possible  with  noninvasive  ultrasound. 
B-mode  ultrasound  imaging  presents  a  display  with  one  axis  representing 
azimuth  and  the  other  axis  representing  range  (Institute  of  Electrical  and 
Electronic  Engineers,  Inc.,  1982,  p.  10).  The  transducer  is  placed  against  the 


body  with  an  impedance  matching  gel  to  insure  good  transmission  between 
the  probe  and  a  vessel.  This  technique  also  can  image  aneurysm  and  thrombi, 
which  may  not  always  be  possible  with  x-ray  (Wells,  1982,  p.  178).  There 
axe  disadvantages  associated  with  this  technique  also.  The  vessel  is  often  not 
easily  imaged  because  of  intervening  acoustically  opaque  tissue,  such  as  bone 
or  the  air  in  lungs.  Additionally,  attenuation  becomes  significant  for  imaging 
vessels  farther  away  from  the  surface,  especially  for  higher  frequency  probes 
that  allow  better  resolution  images.  The  images  of  atherosclerotic  lesions  by 
B-mode  scan  have  had  limited  success  thus  fax  in  characterizing  biological 
composition  and  surface  of  a  lesion. 

Intravascular  techniques  are  a  more  recent  development.  As  well  as  usage 
in  diagnosis  they  are  often  used  in  conjunction  with  atherectomy,  balloon 
angioplasty  and  laser  angioplasty.  There  are  two  intravascular  techniques 
currently  being  developed:  fiberoptic  angioscopy  and  intravascular  ultrasound 
(Coy  et  al.,  1991,  pp.  1811-1812).  Both  of  these  techniques  are  delivered  by 
catheter  through  the  arteries. 

Fiberoptic  angioscopy  provides  visual  evidence  of  thrombus  in  acute  is¬ 
chemia  syndromes,  areas  of  blood  deficiency.  It  also  allows  the  inspection  of 
the  vessel  for  grafting  suitability  and  placement.  Additionally,  the  surface 
of  the  intima  can  be  inspected  and  any  intraluminal  flaps  can  be  detected. 
Unfortunately,  these  imaging  capabilities  can  only  be  done  while  flushing  the 
area  of  any  blood  with  clear  saline.  Another  disadvantage  is  the  inability  to 
image  beyond  the  intima. 

Intravascular  ultrasound  uses  two  kinds  of  transducer  arrangements: 
phased  array  and  rotating  mirror.  Each  system  views  the  structure  of  an 


6 


artery  in  cross-section  including  intima,  media,  and  adventitia,  as  well  as 
some  external  tissue. 

The  images  of  these  layers  bear  a  relationship  to  the  structure  of  the 
artery.  The  image  of  a  normal  artery  shows  a  bright  return  from  the  intima 
and  the  internal  elastic  lamina.  The  next  layer  is  the  non-echogenic  media, 
which  appears  as  a  dark  region.  A  third  layer,  bright  again,  corresponds  to 
the  adventitia  (Yock  et  al.,  1991  p  39B).  This  approach  appears  to  produce 
best  images  in  muscular  arteries,  which  are  generally  the  ones  most  often 
found  with  atherosclerosis.  The  thicknesses  of  these  regions  have  shown  good 
agreement  with  histologic  measurement  with  the  exception  of  the  single  in¬ 
tima  cell  layer  because  it  is  smaller  than  the  resolution  of  most  intravascular 
ultrasonic  imaging  catheters. 

Because  of  the  close  proximity  to  the  lesions  very  high  frequencies  (20  -  30 
megahertz)  can  be  used  with  the  attendant  higher  resolution  without  loss  due 
to  attenuation.  This  also  overcomes  the  disadvantage  of  noninvasive  ultra¬ 
sound  of  intervening  air  and  bone  structures.  By  moving  the  catheter  along 
the  artery,  three-dimensional  structure  can  be  seen.  This  allows  the  presence 
and  structure  of  atherosclerosis  to  be  defined  both  beneath  the  surface  and 
along  the  vessel.  The  limitation  of  the  diameter  of  the  vessel  is  present,  but 
catheters  as  small  as  1.16  millimeter  are  used  (Coy  et  al.,  1991,  p.  1812).  This 
limitation  restricts  the  use  to  lumen  smaller  than  this  diameter  and  means 
that  the  normal  diameter  of  the  artery  will  have  to  be  greater  depending  on 
the  degree  of  stenosis.  However,  because  of  the  high  frequencies  (20  to  30 
megahertz),  imaging  can  take  place  with  the  catheter  nearly  touching  the 
intima.  The  anticipated  use  of  intravascular  ultrasound  is  for  diagnosis  as 
well  as  positioning  of  angioplasty  devices,  and  assessment  of  restenosis  after 
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their  use.  While  there  are  many  advantages  there  are  also  some  disadvan¬ 
tages.  The  scanning  is  perpendicular  to  the  axis  of  the  vessel  which  prevents 
a  view  of  what  possible  obstructions  might  be  ahead,  possibly  leading  to  some 
mechanical  trauma  to  the  vessel. 


Purpose  of  Thesis 


In  a  recent  U.  S.  Department  of  Health  and  Human  Services  report  on 
atherosclerosis  a  recommendation  was  made  to  “support  a  critical  appraisal  of 
newer  techniques”  the  first  of  which  is  intravascular  ultrasound  (1991,  p.  127). 
Among  the  goals  set  forth  was  “to  image  the  structure  and  composition  of 
the  atheroma.”  It  is  the  purpose  of  this  thesis  to  investigate  the  possibility  of 
using  ultrasound  backscatter  for  the  characterization  of  atherosclerosis.  The 
assumption  is  that,  as  with  radar  mapping  of  the  ground  from  aircraft,  oblique 
incidence  waves  can  characterize  the  surface  without  being  overwhelmed  by 
pecular  return. 

Externally  this  could  be  done  with  basically  conventional  imaging  ap¬ 
proaches,  the  only  difference  being  an  off  axis  transducer  to  transmit  the 
signal.  Internal  usage  would  require  the  use  of  a  catheter-mounted  array 
or  similar  device.  It  is  proposed  that  a  forward  scanning  transducer  could 
provide  this  capability  intravascularly. 

The  rationale  behind  this  thesis  is  that  varying  degrees  of  roughness  will 
scatter  incoming  energy  over  differing  angular  spread  and  at  different  levels.  A 
similar  principle  has  been  used  in  radar,  sonar,  and  other  mapping  methods. 
Specular  scattering  or  reflection  is  the  most  direct  return  of  energy  from  a 
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surface  or  object.  If  an  acoustic  beam  is  directed  to  a  smooth  surface  perpen¬ 
dicular  to  it,  the  energy  will  return  along  the  same  path.  If  the  surface  is  not 
perpendicular  to  the  beam  the  energy  will  return  along  an  angle  that  is  equal 
to  the  incident  angle.  A  perfectly  smooth  surface  will  scatter  or  reflect  only 
in  the  specular  direction.  As  roughness  increases  scattering  begins  to  spread 
angularly.  It  is  this  energy  that  is  the  basis  for  sonar  and  radar  mapping.  In 
medical  ultrasound  the  specular  return  is  used  for  characterizing  boundary  of 
tissues  of  different  properties.  The  signal  level  returned  from  an  interface  be¬ 
tween  two  inhomogeneous  materials  is  so  high  that  it  is  sometimes  difficult  to 
retrieve  any  characteristic  of  the  boundary.  With  an  off-axis  method  in  which 
the  effect  of  the  main  clutter  is  minimal,  differences  may  become  apparent. 
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Chapter  2 


BACKSCATTER  FROM  ROUGH  SURFACES 


An  analysis  of  scattering  of  various  forms  of  waves  from  surfaces  is  impor¬ 
tant  in  many  fields.  The  most  common  field  is  radar  where  the  backscatter 
of  radar  energy  allows  the  mapping  of  surface  features  for  navigation  and 
mapping,  detection  of  moisture  for  weather  warning  and  prediction,  and  the 
detection  of  aircraft  for  military  as  well  as  commercial  purpose.  Scattering 
is  also  used  as  the  basis  of  one  form  of  X-ray  imaging  in  medicine  (Shung 
et  al.,  1992,  p.  62).  Acoustic  energy  is  used  in  the  oceans  for  detection  of 
various  submerged  objects  from  enemy  submarines  to  schools  of  fish.  Acous¬ 
tic  energy  is  also  used  in  seismic  exploration  for  oil  and  other  studies  of  the 
morphology  of  the  earth.  A  discussion  of  the  theoretical  basis  for  scattered  en¬ 
ergy  from  surfaces  is  necessary  to  understand  whether  atherosclerotic  plaque 
could  be  imaged.  Since  the  plaque  surface  roughness  characteristic  is  virtually 
unknown,  it  has  been  assumed  to  be  a  Gaussian  distribution. 


Theoretical  Background 


Any  discussion  of  scattering  from  rough  surfaces  must  mention  that  re¬ 
search  began  with  Lord  Rayleigh’s  investigation  of  normal  incident  scattering 
from  a  sinusoidal  surface  in  1877  (Volume  2,  1945,  p.89-96).  From  that  be¬ 
ginning  to  more  recent  developments  the  direction  taken  has  involved  a  move 


from  close  formed  solutions  to  those  that  satisfy  random  roughness  condi¬ 
tions.  Ogilvy’s  recently  published  book  has  the  clearest  presentation  of  the 
necessary  theoretical  development  and  will  be  used  as  the  primary  reference 
to  what  follows  (1991). 

One  of  the  outcomes  of  Lord  Rayleigh’s  work  was  the  generally  used  term, 
Rayleigh  criterion,  as  a  description  of  a  surface’s  smoothness  or  roughness. 
The  geometry  of  plane  wave  incidence  is  illustrated  in  Figure  2.1.  The  incident 
energy  arrives  at  an  angle,  6\,  from  the  normal  to  the  average  height  plane. 
The  incident  energy  that  scatters  in  the  same  direction,  62,  from  two  different 
points  will  have  some  change  in  phase,  A <f>. 

A <f>  **  k[(h\  —  /»2)(co8  0i  +  cos  #2)  +  (X2  —  xi)(sin0i  —  sin #2)]  (2.1) 

The  term,  k,  is  the  wave  number  (2x/A;  A,  wavelength)  of  the  incident  wave. 
If  A<j>  is  very  small  compared  to  *,  the  scattered  waves  will  be  nearly  in 
phase  and  will  constructively  interfere.  If  A<t>  is  approximately  equal  to  x, 
the  scattered  waves  will  be  out  of  phase  and  will  destructively  interfere.  The 
Rayleigh  criterion  limits  A<f>  to  values  less  than  x/2  for  smooth  surfaces.  A 
value  larger  than  x/2  indicates  a  rough  surface  (Ogilvy,  1991,  p.  3).  Replacing 
the  difference  in  h  with  the  average  root  mean  square  (RMS)  surface  deviation, 
a,  the  criterion  becomes: 

Ra  =  kacosdi  <  */4.  (2.2) 

The  symbol,  Ra,  is  called  the  Rayleigh  parameter.  In  the  general  literature 
for  material  roughness  this  same  symbol  is  used  for  the  average  roughness  of 
a  surface,  usually  a  machined  surface  (Bennett  and  Mattsson,  1989,  p.  39). 
In  this  thesis  average  roughness  will  be  symbolized  by  Ra- 


Figure  2.1.  Rough  scattering  from  two  points. 


In  the  discussion  of  roughness  an  important  concept  is  that  roughness  is 
not  an  intrinsic  characteristic  of  the  surface.  A  surface  can  be  considered 
smooth  or  rough  depending  on  the  wavelength  of  the  propagating  energy  in 
the  fluid  medium.  This  can  be  seen  in  the  mutual  relationships  of  the  three 
variables  in  equation  2.2.  The  combined  value  of  ko  cos  &i  determines  whether 
a  surface  is  considered  rough.  In  cases  where  ka  cos  0\  is  nearly  jt/4  the  value 
of  9  may  be  able  to  shift  the  Rayleigh  parameter  from  a  rough  designation  at 
normal  incidence  to  a  smooth  designation  at  nearly  grazing  incidence  (Ogilvy, 
1991,  p.  4). 

Another  area  of  concern  is  the  phase  relationship  of  scattered  wave  fronts 
from  an  extended  surface.  From  equation  2.1  for  a  smooth  surface  (h\  =  hi), 
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the  phase  varies  over  different  values  of  (x2  —  *i).  For  the  case  of  a  smooth 
surface,  the  result  is  the  familiar  beam  pattern  that  clusters  most  of  the 
energy  in  the  specular  direction,  where  0\  =  62-  Additionally,  all  the  energy 
is  in-phase,  constructively  interfering.  As  a  surface  becomes  less  smooth  this 
relationship  becomes  less  clearly  aligned  in  the  specular  direction  and  the 
beam  pattern  spreads  in  angle  and  becomes  less  powerful  in  the  specular 
direction.  This  reduction  is  approximated  by  e-*/2  where: 

g  =  4 k2<r2  cos2  =  4 R2  (2.3) 

Additionally,  as  noted  in  equation  2.1,  the  scattering  in  non-specular  di¬ 
rections  has  some  phase  relationship  that  is  determined  by  the  variations 
across  the  surface.  Due  to  the  random  nature  of  the  surface  fluctuations  the 
phase  relationship  across  a  surface  will  be  distributed  across  all  possible  val¬ 
ues.  The  total  pressure  level  from  this  scattering  will  then  have  an  integrated 
amplitude  that  sums  to  zero.  However,  at  any  angle  there  will  be  some  energy 
scattered  which  has  a  phase  that  is  out  of  phase  with  the  coherent  field  of 
the  specular  scattering.  This  energy  is  termed  the  diffuse  or  incoherent  field 
(Ogilvy,  1991,  pp.  4-6). 

Rough  surfaces  by  their  very  nature  can  not  be  represented  by  any  closed- 
form  description.  This  requires  that  the  surfaces  be  described  by  statistical 
properties.  To  be  able  to  use  statistical  theory  in  the  analysis  of  scattering 
from  rough  surfaces  several  conditions  must  be  satisfied. 

It  is  often  assumed  the  surface  is  isotropic.  This  condition  says  that 
the  surface  characteristics  are  the  same  in  any  direction  chosen.  While  this 
condition  is  not  essential,  it  does  allow  some  simplification  of  the  mathematics. 
A  necessary  condition  for  the  use  of  statistical  theory  is  that  the  surface 
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exhibit  stationarity.  This  condition  is  that  the  statistics  relating  to  one  area 
are  identical  to  any  other  area  of  the  same  surface.  The  probability  that  a 
surface  has  an  expected  value  for  height  at  one  point  does  not  change  as  the 
point  is  changed.  The  final  condition  for  use  of  statistical  theory  is  ergodicity. 
This  is  the  condition  that  any  stationary  statistical  surface  will  retain  the 
same  properties  regardless  of  looking  at  a  single  sample  or  an  ensemble  of 
samples  (Bendat  and  Piersol,  1971,  pp.  10-12). 

There  are  two  statistical  values  that  can  define  a  surface.  The  first  is 
the  variation  of  heights  along  a  surface  or  the  surface  height  distribution. 
The  second  value  is  the  spread  of  heights  along  the  surface  or  the  surface 
correlation  function  (Ogilvy,  1991,  pp.  9-10). 

Suppose  that  the  surface  height  is  h(r),  where  h  is  the  height  above  or 
below  a  reference  surface,  and  r  is  the  position  vector  on  the  surface.  For  this 
random  surface  there  is  a  statistical  height  distribution,  p{h),  where  p(h)dh 
is  the  probability  that  a  point  on  the  the  surface  has  a  height  between  h  and 
k  +  dh  away  from  the  mean  surface.  The  mean  surface  is  usually  defined  as: 

{h)t  *  [°°  hp(h)dh  =  0  (2.4) 

J—OO 

where  ( h)a  represents  the  expected  value  or  the  average  across  a  surface.  The 
root  mean  square  height  is  the  standard  deviation  and  is  defined  as 

a  =  y/(h?)l.  (2.5) 

Surface  roughness  can  also  be  calculated  by  the  arithmetic  mean  or  center 
line  average  (cla)  which  is  defined  as 

Rda  =  /  \h\p(h)dh. 

J—OO 


(2.6) 


Considerable  argument  about  the  aptness  of  a  Gaussian  assumption  has 
been  made  (Ogilvy,  1991,  pp.  11-12).  In  general,  the  conclusion  is  that  it  is 
a  good  assumption  until  proven  otherwise.  For  a  Gaussian  distribution  the 
formula  is: 

For  a  given  Gaussian  distribution  p(h),  the  center  line  average  has  the  follow¬ 
ing  relationship  to  the  standard  deviation  (Ogilvy,  1991,  p.  10): 

Rda  =  ey/V*.  (2.7) 

Measurements  of  surface  roughness  can  be  done  by  several  methods.  The 
usual  method  is  average  roughness,  Ra,  and  is  done  with  a  profilometry  de¬ 
vice.  This  device  moves  a  stylus  across  the  surface  and  measures  the  deflec¬ 
tions  as  it  moves.  An  average  surface  roughness  is  calculated  and  the  distance, 
hi,  above  or  below  this  surface  is  summed  and  then  divided  by  the  length  L 
(Talysurf  10  Operator’s  Manual,  1977,  p.  3). 

RA=^^__y:^  (28) 

The  definition  for  RA  is  identical  to  the  center  line  average.  The  surface 
roughness  is  also  described  as  the  root  mean  square  roughness  a.  This  is 
defined  as  the  square  root  of  the  mean  value  of  squared  heights  from  the 
mean  surface  level  (J.  M.  Bennett  and  L.  Mattsson,  1989,  p.  38). 


The  values  for  Ra  and  a  are  similar  if  there  are  no  large  deviations  from  the 
mean  surface  level.  If  there  are  many  large  deviations  they  will  dominate  the 


sum  and  a  will  be  larger  than  R&  (J.  M.  Bennett  and  L.  Mattsson,  1989,  p. 
39).  While  Ogilvy  says  there  is  a  scale  factor  of  yj2 /n,  the  values  calculated 
or  measured  for  R a  are  assumed  equivalent  to  a  in  this  thesis. 

The  second  value  to  define  the  rough  surface  is  correlation  length.  The 
RMS  height  does  not  distinguish  between  a  rapidly  changing  surface  and  a 
slowly  varying  one.  To  find  this  length  scale  factor  a  correlation  function  must 
be  used  (Ogilvy,  1991,  p.  13).  The  autocorrelation  function  is  defined  as: 


D{R)  =  (h(r)h(r  +  R))a 


(2.10) 


Tb?s  is  the  normalized  version  of  the  autocovariance  function  which  lacks  the 
<r2  term.  The  autocovariance  function  in  summation  notation  is: 

i  N-l 

m  =  T7  E  Mi+i.  /  *  0, 1, 2, . . . ,  N  -  1  (2.11) 

iV  i=l 

The  values  of  A,  are  measured  from  the  mean  surface  (Bennett  and  Mattsson, 
1989,  p.  46).  If  a  surface  is  Gaussian  the  correlation  function  is  given  by 

D(R)  =  exp(-^).  (2.12) 

Here  the  factor  A0  is  the  correlation  length.  The  correlation  length  is  defined 
as  the  distance  along  the  correlation  function  where  the  value  drops  to  e-1 
times  the  value  at  zero.  Notice  that  the  value  R  in  equation  2.12  is  no  longer  a 
vector.  This  is  because  the  surface  is  assumed  to  be  isotropic.  The  direction 
of  correlation  measurement  will  yield  the  same  result  for  isotropic  surfaces 
regardless  of  direction  chosen. 

An  alternative  description  of  the  surface  is  to  look  at  the  power  spectrum. 
This  is  the  Fourier  transform  of  the  un-normalized  correlation  function. 


p<k)  =  (^/->(R)ei(lMRM(R) 


(2.13) 


For  a  Gaussian  surface  with  am  anisotropic  surface  the  power  spectrum  has 
the  following  form 


P(k  i,k2)  = 


?lhhexp{JM)exp{JM 

An  A  4 


(2.14) 


The  spectrum  is  also  Gaussian  with  a  mean  of  zero  and  a  standard  deviation 
of  y/2/\j  where  j  takes  the  values  1,2  for  the  x  and  y  directions  (Ogilvy, 
1991,  pp.  17-18).  If  the  surface  is  isotropic  the  formula  simplifies  to: 


P(k)  = 


<T2  \l 
An 


(2.15) 


There  are  some  additional  functions  that  are  used  in  the  description  of 
surface  roughness.  While  they  are  not  used  in  this  thesis  they  are  mentioned 
for  completeness.  One  such  function  is  the  surface  structure  function.  This 
function  is  the  mean  square  of  the  difference  in  height.  This  function  has  the 
advantage  of  being  independent  of  the  reference  height  chosen.  Unfortunately 
it  does  not  figure  directly  in  the  theory  of  wave  scattering  (Ogilvy,  1991, 
p.  17).  Another  function  is  the  characteristic  function  which  is  the  Fourier 
transform  of  the  probability  density  function. 


xW  =  f°°  p(h)'Mdh 

J—oo 


(2.16) 


This  has  the  advantage  of  including  the  phase  modulation  of  the  wave  at 
a  rough  surface.  In  addition  to  these  two  functions  there  are  higher  order 
surface  functions.  The  first  is  the  two  point  height  probability  distribution.  If 
two  sets  of  random  variables  are  independent  their  joint  probability  becomes 
the  product  of  their  single  point  probability  functions.  There  are  also  higher 
order  surface  correlations.  Two  such  functions  are  skewness  and  kurtosis. 
Skewness  is  a  third  moment  measure  of  how  the  random  surface  has  more 
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peaks  or  valleys.  A  positive  skewness  indicates  that  the  average  surface  level 
has  been  calculated  higher  than  it  should  be.  Kurtosis  is  a  fourth  order 
moment  description  of  the  surface.  A  kurtosis  value  of  3  indicates  a  Gaussian 
distribution  of  height.  Less  than  three  indicates  a  surface  with  fewer  extremes 
than  a  Gaussian  surface.  A  kurtosis  greater  than  3  indicates  a  surface  with 
more  small  height  variations  than  a  Gaussian  surface  (Bennett  and  Mattsson, 
1989,  p.  43).  Finally  there  is  the  slope  distribution  function.  This  measures 
the  slope  of  the  surface  rather  than  the  height  (Bennett  and  Mattsson,  1989, 
p.  47). 

There  are  two  approximations  used  in  the  theory  of  scattering  from  rough 
surfaces.  They  are  perturbation  theory  and  tangent  plane  or  Kirchhoff  theory. 
In  general,  they  cover  different  conditions  but  do  have  some  overlap.  For 
both  theories  the  acoustic  field  is  assumed  to  be  composed  of  the  incident 
and  scattered  fields. 

^(r)  =  V>,we(r)  +  ^c(r)  (2.17) 

The  function  ip  is  usually  the  velocity  potential  for  acoustic  waves. 

Perturbation  theory  is  used  for  slightly  rough  surfaces.  There  are  restric¬ 
tions  in  its  use: 

Jb|A(*,y)|<l  (2-18) 

|Vfi(x,y)|  <  1  (2.19) 

These  restrictions  arise  from  the  use  of  Taylor  series  expansion  on  the  mean 
scattering  surface.  Equation  2.18  indicates  that  the  surface  height  must  be 
much  less  than  the  wavelength.  Equation  2.19  indicates  that  the  slopes  are 
very  small,  indicating  a  gently  sloping  undulation  where  phase  relationships 
are  generally  preserved  (Ogilvy,  1991,  pp.  38-39). 
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The  Kirchhoff  theory  makes  an  approximation  to  the  field  on  the  surface 
of  the  scatterer.  Each  point  on  the  surface  is  assumed  to  be  part  of  an  infinite 
plane.  As  Ogilvy  notes  the  theory  is  exact  for  an  infinite  smooth  plane,  but 
becomes  an  approximation  for  scatterers  that  are  finite,  rough  or  non-planar 
(1991,  p.73).  The  choice  of  theory  becomes  one  of  selecting  a  theoretical  base 
that  is  at  best  an  approximation  but  at  least  is  valid  in  the  regions  of  surface 
roughness  relative  to  the  wavelength  of  interest. 


Development  for  Application 


The  Kirchhoff  theory  for  acoustic  waves  is  developed  for  the  application  of 
modelling  surfaces  of  atherosclerotic  plaque.  Using  a  Helmholtz  integral  scat¬ 
tering  formula  to  compute  the  scattered  field  results  in  the  following  formula 
starting  from  equation  2.17. 

V’(r)  =  ~  G(r. (2.20) 


In  equation  2.20  S0  is  the  surface  of  the  scatterer,  r  is  in  a  closed  volume 
containing  no  sources  from  rp,nc(r)  or  0sc(r)  or  outside  a  closed  volume  that 
encloses  all  the  sources.  The  unit  surface  normal,  n0,  points  towards  the 
source.  The  function,  G( r,  rc),  is  the  acoustic  Green’s  function  of  the  field 
at  r  for  a  source  at  r„.  The  Green’s  function  for  full  space  scattering  from  a 
finite  surface  is: 


G(r,r0)  = 


txp{ik\ r  —  r0[) 


(2.21) 


4jr|r  -  r0| 

Here,  r0,  is  on  the  scatterer  surface  and  r  is  some  distance  from  the  surface. 
When  the  surface  S0  is  closed  the  scattered  field,  rfrac(r0) ,  in  equation  2.20 
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is  interchangeable  with  total  field,  xp( rc).  This  is  due  to  no  sources  of  ip'nc 
within  S0. 


/  [0,nc(ro) 

•'So 


o  na 


dnB 


(2.22) 


The  choice  of  surface  depends  on  the  field  to  be  used  in  the  integrand  of 
equation  2.20.  If  the  scattered  field  is  used  the  rough  surface  must  connect 
with  a  surface  at  infinity.  The  surface  at  infinity  has  a  scattered  field  of  zero. 
If  the  total  field  is  used  in  equation  2.20,  the  surface  is  closed  with  a  surface 
just  below  the  rough  surface  with  a  total  field  of  zero  on  this  closing  portion. 
The  first  assumption  ignores  the  edges  of  finite  surfaces.  The  second  surface 
integral  will  include  the  effect  of  diffraction  from  the  edges.  Ogilvy  citing 
Beckmann  and  Spizzichino,  using  the  second  method  (  otal  field)  gives  the 
scattered  field  as  (1991,  p.77): 


*"(r)  =  *(r)  -  0  (r) 


=  f  Mro) 


dG(r,T0) 

dna 


9rp(r0), 

dna 


—  G(r,ro)—^~]dS0 


(2.23) 


Using  the  reflection  coefficient  for  plane  waves,  R0,  the  total  field  on  the 
surface  is  (the  Kirchhoff  approximation): 


V>(r0)  =  [1  +  #o(r0)]^mc(r0) 


Assuming  a  plane  monochromatic  wave  of  the  form 


(2.24) 


0,nc(r)  =  e‘k,"'r 


(2.25) 


where  ki„c  is  the  incident  wavevector.  This  sets  the  normal  derivative  as 


dtHr q) 
dna 


=  i[l  -  i2o(ro)]ktoc  •  riptnc(r). 


(2.26) 
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The  next  step  is  to  make  the  far  field  approximation,  r  >  r0.  Additional 
requirements  are  often  placed  on  the  problem  for  incident  beams:  r  d2 /A0 
and  d2  A2  where  d  is  the  diameter  of  the  illuminated  surface.  This  latter 
restriction  is  necessary  also  for  the  stochastic  problem.  The  far  field  approx¬ 
imation  results  in  the  scattered  field: 

tT(r)  =  ^  /  [(*.k-  -  k+)  •  no]e,k"r“dS0  (2.27) 

47rr  JSo 

where  the  wavevectors  are  defined  as: 

k  =  kjae  kae 


k+  =  kinc  +  k^ 


(2.28) 


The  integration  is  over  the  surface,  S0,  which  undulates.  Converting  the 
integration  over  the  mean  plane  surface  requires  that  the  element  area  become: 


nodS0  ~  -  j T7—  +  &)dSM 

OXq  oyQ 


(2.29) 


since  the  unit  normal  is 


n.  =  +  t  (2.30) 

yf{l+(dhldX'f  +  (dhldy.f) 

This  change  also  requires  the  additional  assumption  that  no  point  on  the 
surface  has  infinite  gradient  (Ogilvy,  1991,  p.79). 

For  the  standard  geometry  shown  in  Figure  2.2,  where  6\  is  the  angle  from 
the  z  axis  to  the  incident  wave  vector  which  projects  onto  the  x  axis,  62  is  the 
angle  from  the  z  axis  to  the  scattered  wave  vector,  and  63  is  the  angle  from 
the  x  axis  to  the  projection  of  the  scattered  wave  vector  onto  the  xy  plane, 
the  incident  and  scattered  wave  vectors  become 


kinc  =  fc(i  sin  $i  —  k  cos  61) 


kte  =  &(i  sin  &2  cos  $3  ■+■  j  sin  62  sin  #3  4-  k  cos  #2) 


The  scattered  field  becomes 


ieikr  r 

^*c(r)  =  — —  /  ( adh/dx0  +  bdh/dy0  —  c) 

4xr  JSj£ 


xexp[ik(Ax0  +  By0  +  Ch(x0,  yo)]dx0dy0. 


The  respective  variables  are  as  follows: 


A 

A  =  k~  =  sin  0\  —  sin  62  cos  63 

B  =  k~  =  —  sin  02  sin  63 
C  —  kj  =  —  (cos  +  cos  62) 

a  =  (k+  —  Rok~)x  —  sin  (1  —  Ro)  -f  sin  $2  cos  03(1  -f  Ro) 
b  =  (k+  —  Rok~)t  =  sin  02  sin  £3(1  —  Ro) 
c  —  (k+  —  Rok~)t  =  cos  02(1  +  Ro)  —  cos  0i  (1  —  Ro) 

An  additional  step  can  remove  the  gradient  expressions  by  the 
integration  by  parts  and  assuming  the  area  of  integration  extends  over 
x0  <  X  and  —Y  <  y0  <  Y .  The  scattering  function  now  has  the  form: 

*“(r)  =  +  fr  +  C)  JSu  +  *, 

where  the  edge  effects,  Vv, 

rl>t e  =  —7— -[7-;  [(eik*x*°)  -  eik^~x^)dy0 
4 at  fc(7  ./ 
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(2.31) 


(2.32) 


(2.33) 

use  of 
-A  < 

(2.34) 


y  (c‘**(*o,r)  _  e.fc0(*o,-y'))<f;ro]j 


(2.35) 
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Figure  2.2.  Coordinate  system. 

and  the  phase  is 

y0)  =  Ax0  +  By0  +  Ch(x0 ,  y0)-  (2.36) 

The  angular  terms  can  be  consolidated  into  one  function 

m, »»■»«)  =  !(£  +  f?  +  c),  (2.37) 

yielding  (Ogilvy,  1991,  p.  83) 

<T(r)  =  ^-2F(«i,  «2,  «s)  J  (2.38) 

The  edge  effects  are  often  assumed  negligible.  While  this  is  a  convenient 
assumption  it  is  only  true  for  the  near  specular  case  (0i  =  #2,  #3  =  0)  (Ogilvy, 
1991,  p.  85). 
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This  formula  is  the  general  form  and  now  can  be  used  to  look  at  the 
coherent  and  diffuse  fields.  With  the  introduction  of  the  probability  density 
function  the  average  mean  acoustic  field  becomes 

(V>'c)  =  Z1S—2F(0u$2,03)  f  [°°  eik^Xo^p(h)dhdx0dy0  +  <0e>  (2.39) 
4jt  r  JS J—oo 

The  amplitude  of  the  coherent  field,  after  some  manipulation,  becomes 

W“)  =  =^.AucxikC)(^^)(S-^^)  =  X(*C»1‘  (2.40) 

where  tpgC  is  the  field  scattered  from  a  smooth  surface  of  area  A\f.  For  a 
Gaussian  distribution  the  characteristic  function  is 

(2-41) 

From  this  equation  2.40  becomes 

<<r> = x(*c)0r = (2.42) 

where 

g  =  jfc2C2<r2  (2.43) 

For  specular  reflection  the  Rayleigh  parameter,  R+,  is  equivalent  to  \y/g- 
For  the  diffuse  field  the  mean  amplitude  would  sum  to  zero  but  the  average 
intensity,  (/<*),  can  be  evaluated. 

(Id)  =  (P't")  -  «■“><?“)  (2-44) 

Three  additional  assumptions  are  made  about  the  surface.  The  first  is  that  the 
dimensions  of  the  surface  are  larger  than  the  correlation  length.  The  second 
is  that  the  surface  of  interest  is  isotropic  in  order  to  use  spherical  coordinates. 
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The  third  is  that  the  edge  effects,  (t/>e),  are  not  stochastic  and  give  no  con¬ 
tribution  to  the  diffuse  field  (A(±A,±Y'’)  =  0)).  From  these  assumptions  the 
following  is  true  (Ogilvy,  1991,  p.  87) 

Jfc2  p2  ,00  / - 

{Id)  =  l  Jo^x2  +  y2) 

X[X2 (kC,  - kC ,  R )  -  x(kC)^kC)]RdR.  (2.45) 

From  Ogilvy  the  two  dimensional  characteristic  function  is  (1991,  p.  89) 

X2 (kC,  -kC,  R)  =  exp(-(k2C2<72[l  -  D(R)}))  (2.46) 

and  the  surface  correlation  function  D(R)  is 

D(R)  =  e-R2/x*.  (2.47) 


Finally  this  yields  an  expression  for  the  diffuse  field  intensity  of 


k2F2  \2c~’ 


o°  gn 


k2(A2  +  B2)X2 


4n 


) 


(2.48) 


This  carries  the  restriction  that  g  <  1  for  the  series  to  converge.  The  total 
field  intensity  is  the  sum  of  the  coherent  and  diffuse  fields. 


(/)  =  he-’  +  (h)  (2.49) 

The  result  of  equation  2.49  is  that  fields  of  scattered  energy  from  various 
Gaussian  roughness  can  be  predicted  if  the  RMS  height,  wavevector,  incident 
and  scattering  angles,  illuminated  area  and  correlation  length  are  known. 
Ogilvy  divides  surfaces  into  three  regimes  for  evaluation  (1991,  pp  91-93).  The 
divisions  are  g  1,  for  slightly  rough  surfaces,  g  ~  1,  for  moderately  rough 
surfaces,  and  g  >  1,  for  very  rough  surfaces.  The  first  and  second  regions 
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give  information  about  the  diffuse  fields.  The  third  region  is  dominated  by  the 
diffuse  field  and  it  becomes  the  total  field.  Additionally,  in  the  second  region 
the  diffuse  field  is  described  as  a  range  between  a  minimum  and  maximum 
possible. 

For  the  slightly  rough  surface  the  diffuse  field  intensity  is  the  first  term 
in  equation  2.48. 


,,,  ...  ,  t2(A2  +  fi2)^ 

(li)  ~  V  'A«eip( - ; - ) 


4xr2 


kAC2  F 2 

-±J-AMP(kA,kB) 


(2.50) 


where  P(kA,  kB)  is  equation  2.14. 

For  the  moderately  rough  surface  the  range  of  diffuse  intensity  is  in  the 
range: 


k2F2  A2 
4xr2 


u 2  f2x2 

~^AU  (2.51) 


For  very  rough  surfaces  the  total  intensity  magnitude  has  the  following 
form: 

A  B 

{!)  =  -2  "-AjfPntp,  (2.52) 


where 

,  ^  /  *2(*2  +  y2h 

Pi2(ar,y)  =  4^2cx*K - £5 - )  <2-53) 

From  equations  2.51,  2.52,  and  2.53  some  projections  can  be  made  about 
what  level  of  backscatter  will  be  seen  given  an  RMS  roughness  and  a  correla¬ 
tion  length. 


26 


Chapter  3 


MATERIALS 


In  order  to  study  the  feasibility  of  imaging  atherosclerotic  plaque  and 
quantitating  its  surface  roughness  suitable  samples  are  needed  that  duplicate 
the  surfaces  seen  in  vivo.  While  real  tissue  samples  can  be  utilized,  it  is  more 
convenient,  as  well  as  more  accurate,  to  use  samples  of  a  known  roughness  and 
correlation  length.  The  references  found  in  this  research  contain  no  data  of 
RMS  surface  roughness  or  correlation  length  for  typical  lesions.  Because  data 
on  the  roughness  of  atherosclerotic  plaque  is  only  of  an  anecdotal  nature,  a 
wide  range  of  roughness  samples  is  required.  The  samples  of  roughness  are  in 
three  forms.  The  first  samples  are  Plexiglas  blocks  that  have  been  sanded  with 
various  grades  of  sandpaper.  The  second  samples  are  also  Plexiglas  blocks 
with  glass  beads  or  carbon  particles  bonded  to  their  surfaces.  The  third  type 
is  a  cast  surface  comparator  from  the  Aluminum  Association  (300  19th  Street, 
N.W.,  Washington,  D.C.).  In  addition  to  the  materials  mentioned,  there  are 
three  tissue  samples  used  for  demonstration  purposes. 

These  samples  provide  an  opportunity  to  test  the  scattering  theory  with 
real  surfaces.  There  are  some  dissimilarities  between  the  sample  surfaces  and 
blood  vessel  tissue  however.  The  basic  difference  is  that  the  materials  have 
different  density  and  sound  speed  than  arterial  tissue.  This  is  discussed  at 
the  end  of  this  chapter. 
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Sanded  Surfaces 


The  Plexiglas  blocks  are  5  centimeter  square  surfaces,  2  centimeters  thick. 
They  are  illustrated  in  Figure  3.1.  The  holes  through  the  sides  are  for  mount¬ 
ing.  The  Plexiglas  blocks  are  sanded  with  the  following  sandpaper  types:  60, 
100,  180,  280,  and  400  grit  roughness.  The  sanding  is  accomplished  by  apply¬ 
ing  the  paper  to  the  surface  in  the  same  direction  until  the  gloss  of  the  smooth 
surface  is  replaced  with  a  uniform  sheen.  This  sheen  is  an  indication  that  the 
abrasive  has  removed  as  much  of  the  Plexiglas  as  is  possible  but  before  it  has 
begun  to  actually  remove  the  highest  peaks. 


Figure  3.1.  Plexiglas  block. 
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Sandpaper  is  usually  described  by  a  grit  number.  While  often  the  grit  is 
assumed  to  represent  the  number  of  particles  per  linear  inch,  this  is  not  the 
case.  This  number  is  a  reference  to  the  size  of  the  sieve  used  to  separate  the 
abrasive  particles.  The  grit  number  specifies  the  number  of  holes  per  linear 
inch  in  the  sieve.  The  following  Table  3.1  shows  the  actual  sizes  of  the  grains, 
and  the  associated  grit  number  for  closed  coat  sandpaper  (Norton  Co.,  1990, 
p.  3).  (Closed  coat  refers  to  the  fact  that  the  grains  completely  cover  the 
backing  paper.  Open  coat  sandpaper  has  only  50%  to  60%  of  the  backing 
<x>vered.)  As  can  be  seen  in  the  table  the  sanded  surface  roughness  is  much 
less  than  the  particle  size  of  the  sandpaper  used  to  makes  it.  This  is  due  to 
the  distribution  of  the  particles  on  the  sandpaper.  Adjacent  particles  prevent 
the  full  penetration  into  the  surface  of  any  given  particle. 


Table  3.1.  Sandpaper  information. 


Grit 

Number 

Particle 

Size 

pm 

Ra 

ftm 

fxm 

268.0  ^ 

15*  "I 

~~  lTr* 

■■miSSJBHH 

9.5 

180 

Y71 

7.5 

IHHaSflH 

jUgpEgJMMB! 

23.6 

4.1* 

*  calculated 


Three  of  the  samples  (100,  180,  and  280  grit)  were  tested  on  a  Talysurf  10 
profilometer.  The  profilometer  is  only  capable  of  measuring  small  roughness 
values.  Anything  above  a  value  of  ~  3.0 fxm  would  give  inaccurate  results 
and  perhaps  damage  the  machine.  The  profilometer  moves  a  stylus  over  the 
surface  and  measures  the  height  of  the  surface  relative  to  a  mean  surface. 
The  machine  also  makes  a  trace  of  the  surface  on  a  strip  chart.  A  linear 
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regression  for  the  roughness  measured  on  the  profilometer  to  the  grain  size 
of  the  sandpaper  was  made.  The  correlation  for  the  regression  is  0.999.  The 
formula  for  the  linear  regression  is: 

Ra  =  0.30137  +  0.01 196d,  (3.1) 


where  d  is  particle  size. 

The  correlation  length  for  the  sanded  surfaces  is  found  from  the  Talysurf 
10  strip  charts.  Sixty-four  points  are  sampled  at  5 nm  from  the  chart  and 
then  correlated  and  the  correlation  length  found.  The  measured  values  are 
then  correlated  with  particle  size,  grit,  and  roughness.  The  correlation  length 
had  a  correlation  with  particle  size  of  —0.986,  with  grit  of  0.993,  and  with 
roughness  of  —0.993.  The  values  predicted  from  the  grit  correlation  are  used 
for  the  two  surfaces  not  measured.  The  formula  for  the  linear  regression  is: 

X0  =  2.609  +  0.025$rt*,  (3.2) 

The  predicted  values  are  marked  with  an  asterisk  in  Table  3.1. 

In  addition  to  the  sanded  surfaces,  a  block  made  with  the  400  grit  paper 
has  five  regions  of  extreme  roughness  of  differing  areas.  The  areas  are  2,  3,  4,  5, 
and  10  millimeter  squares.  This  is  a  resolution  sample  that  is  for  maximum 
contrast.  The  roughened  regions  were  made  with  a  stylus  and  by  physical 
inspection  are  rougher  than  any  of  the  samp1'"  s  available.  Unfortunately  no 
RMS  roughness  or  correlation  length  are  known.  This  resolution  sample  is 
shown  in  Figure  3.2. 
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Glass  Beads 


The  glass  bead  surfaces  are  prepared  by  placing  glue  on  the  portion  of  the 
Plexiglas  for  the  beads.  The  beads  axe  then  spread  over  the  area  in  sufficient 
quantity  to  assure  complete  coverage. 

The  glass  bead  surfaces  required  an  analytical  approach  for  finding  the 
roughness  because  of  the  profilometer’s  sensitivity.  The  calculation  of  the 
roughness  of  a  surface  is  based  on  the  following  formula: 

Ra  =  —  k2  +  *3--—  =t  fL  \h\dL  (3.3) 

L  L  Jo 

The  average  roughness,  Ra,  is  used  generally  for  machined  surfaces  (Bennett 
and  Mattsson,  1989,  p.39).  This  formula  implies  that  a  horizontal  midline 
exists  where  the  area  between  the  line  of  the  surface  trace  above  the  midline 
equals  the  area  of  the  surface  trace  below  the  midline  and  the  midline.  With 
a  few  assumptions  the  roughness  can  be  calculated  for  the  beaded  surfaces. 

The  first  assumption  is  that  the  beads  are  spherical  and  of  the  same 
diameter.  While  there  is  some  variation  among  a  group  of  beads,  they  axe 
reasonably  uniform.  The  second  assumption  is  that  the  beads  on  the  surface 
are  distributed  as  densely  as  possible.  This  is  possible  with  the  uniformity  of 
the  spheres.  The  next  assumption  is  that  the  glue  used  to  apply  the  beads 
to  the  surface  fills  the  volume  below  the  spheres  to  their  equators.  This 
means  the  surface  is  composed  of  hemispheres.  The  next  assumption  is  that 
the  nominal  roughness  can  be  calculated  by  taking  an  average  of  two  trace 
directions.  Because  of  the  dense  packing  of  the  spheres  there  is  a  direction 
where  the  spheres  will  be  touching  at  their  equators.  Also  because  of  the 
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uniformity  of  the  spheres  this  is  a  straight  line.  Because  of  the  symmetry  of 
the  arrangement  this  alignment  will  exist  at  60°  increments  around  a  central 
bead,  as  in  Figure  3.3.  As  the  angle  sweeps  between  the  hexagonal  peaks 
a  maximum  of  roughness  is  found  at  the  midpoint  between  the  60°  spokes. 
A  line  taken  anywhere  else  will  result  in  a  value  between  these  two  values. 
The  average  is  a  nominal  value.  It  is  not  the  exact  value,  but  should  yield  a 
reasonable  value  considering  the  other  assumptions. 


Figure  3.3.  Beaded  surface  from  above. 


The  direction  where  the  spheres  touch  will  be  referred  to  as  A  —  A,  as 
in  Figure  3.4.  The  other  line,  30°  away,  will  be  referred  to  as  B  -  B,  as  in 
Figure  3.5.  Both  lines  are  composed  of  repeating  segments. 
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R 


A  -  A 


Figure  3.4.  Beaded  surface  line  A  —  A. 

For  the  A  —  A  line,  the  segment  length  is  the  same  as  the  diameter  of 
the  sphere,  2 R.  The  roughness  is  calculated  by  summing  the  areas  above 
and  below  a  midline  and  then  dividing  by  the  segment  length.  The  height 
of  the  midline,  h,  is  found  by  equating  the  areas  above  and  below.  The 
area  of  the  semicircle  is  ttR2/2.  The  areas  of  concern,  a,  /?,  and  7,  are 
described  as  the  area  below  h  and  above  the  semicircle,  the  area  above  h  but 
below  the  semicircle,  and  the  area  below  h  after  the  semicircle  drops  below 
h,  respectively.  Areas  a  and  7  are  the  same  size.  To  satisfy  equation  3.3  the 
following  is  true: 

a  +  7  =  /?,  a  =  7  (3.4) 
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B-B 


Figure  3.5.  Beaded  surface  line  B  —  B. 

The  area  of  0  is  the  difference  between  the  total  area  of  the  semicircle  minus 
the  area  below  h  inside  the  semicircle: 

0  =  7TjR2/2  —  (a  +  6  +  c  4-  d),  a  =  b,  c  —  d  (3.5) 

The  figures  show  that  6  describes  the  angle  between  the  base  of  the  semicircle 
and  the  intersection  of  h  and  the  semicircle.  From  this  we  have  the  following 
area  formulas: 

a  =  OR2/ 2,  9  in  radians  (3.6) 

b  =  ^RcosORsinO  (3.7) 

0  =  R2(^r  —9  —  sin  9  cos  9) 


(3.8) 


a  =  7  =  R2(  sin  0  —  6/2  —  cos6sin0/2) 
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(3-9) 

Solving  for  the  angle  9: 

6  =  arcsin  -  (3.10) 

Therefore  for  line  segment  A  —  A  the  roughness,  Ra  is: 

Ra  =  R2(ir/2  +  2  sin  6  -  26  -  2  sin  6  cos  6)  =  0.1813/2  (3.11) 

2/i 

For  the  B  —  B  line  the  repeating  distance  includes  the  space  between  two 
spheres  that  are  oriented  at  a  30°  angle  to  the  A  -  A  line.  This  additional 
area  below  h  is  labelled  6  and  is  summed  with  a  and  7. 


o  +  7  +  *  =  0,  a  =  7  (3.12) 

The  distance  can  be  found  easily  from  geometry.  The  distance  between  the 
centers  of  spheres  on  the  B  -  B  line  is  2y/3R.  Therefore  6  is  defined  by: 

6  =  h(>/3  —  l)R  =  2sin0(%/3  —  l)R2  (3.13) 

From  this  formula  comes  a  new  value  for  6. 


6  =  arcsin 


4 1/5 


(3.14) 


This  gives  a  new  value  for  the  roughness  along  the  B  —  B  line. 


Ra  =  ^j^=R2  (ir /2 +  2y/Z  sin  9  —  26  —  2  sin  9  cos  9)  =  0.4018/2  (3.15) 


Averaging  the  two  values  gives  the  roughness  as  0.2916/2.  Figure  3.6  shows 
a  plot  of  Ra  as  a  function  of  the  spacing  between  beads.  Both  the  roughness 
and  spacing  are  expressed  as  fractions  of  the  radius. 


Ratio  of  Bead  Spacing  to  Radius 


Figure  3.6.  Ra  as  function  of  bead  spacing. 


For  the  sake  of  argument  calculating  the  case  of  no  glue  at  all,  h  is  above 
the  substrate  by  an  increased  distance  R,  6  is  increased  by  2R2(\/3  —  1). 


6  =  2(1+  stn0)(V5  -  l)fl2  (3.16) 

6  =  s'n’1(73(!  +  1  -  <317) 


This  case  yields  then  a  roughness: 

Ra  =  +  2(>/3  -  1)  +  2sin  9 


+2(  'Ji  —  1)  sin  B  —  29  —  2  sin  $  cos  0)  —  0.87 1 3 


(3.18) 


relation  Height  Ratio 


These  results  show  that  a  dense  pack  set  of  spheres  of  uniform  size  will  have  a 
maximum  roughness  associated  with  some  multiple  of  the  radius  that  is  less 
than  one. 


Correlation  length  for  the  bead  surface  can  also  be  calculated.  Figure 
3.7  shows  a  function  depicting  the  bead  surface  and  the  correlation  function. 
There  are  32  points  per  bead  profile  and  the  correlation  drops  below  e_1  at 
the  tenth  noint.  The  correlation  length  is  equal  to  5/16 R. 


Figure  3.7.  Bead  profile  and  correlation  function. 


Two  bead  surfaces  have  been  prepared  and  used.  The  first  surface  has 
the  largest  bead  size  covering  the  entire  surface.  The  surface  distribution 
is  several  layers  deep  due  to  too  much  glue  in  the  application  of  the  beads 


which  allowed  the  beads  to  become  more  than  one  layer  deep.  The  second 
sample  surface  for  beads  is  a  Plexiglas  block  with  5  millimeter  wide  vertical 
strips  composed  of  five  different  sizes  of  spheres.  The  pattern  (237 fim,  50 nm, 
105 fim,  150 fim,  215/um)  is  repeated  twice  for  the  full  width  of  the  block.  The 
data  for  the  beaded  surfaces  are  displayed  in  Table  3.2  and  the  surface  is 
shown  in  Figure  3.8. 


Table  3.2.  Glass  bead  information. 


Dia. 

ftm 

Range 

ftm 

Ao 

fim 

HHKHMHH 

7125 

“  7.81 

15.23 

15741  “ 

2T775 

23.44 

215 

31.18 

33759 

237  ' 

177  -  297 

r~  34.37 

Cast  Surface  Comparator 

The  cast  surface  comparator  is  shown  in  Figure  3.9.  This  surface  is  used 
as  a  standard  for  surface  roughness  in  the  aluminum  casting  industry.  The 
surface  is  divided  into  nine  areas  with  specified  roughness  in  microinches. 
Table  3.3  shows  the  conversion  to  units  of  micrometers.  The  table  also  lists  a 
nominal  correlation  length  that  is  5.4  times  the  Ra  value.  This  is  the  average 
of  the  measured  sanded  surfaces  compared  to  their  Ra  values.  While  this 
is  probably  not  accurate  it  does  provide  a  nominal  value.  Contact  with  the 
Aluminum  Association  Inc.  indicated  that  they  had  no  data  on  the  correlation 
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length  associated  with  their  standard.  The  surfaces  are  far  too  rough  for  the 
Talysurf  profilometer. 


Table  3.3.  Cast  surface  information. 


Ra 

Ra 

fiin 

fim 

fim 

i?5 

L53 

■”  '  8.26  ~ 1 

miiliMi 

1517 

~  3713  “  1 

712 

41.15 

5713 

U722 

7575 

98.77 

122.86 

13314  . 

Tissue  Samples 


Three  samples  of  artery  tissue  have  also  been  studied.  The  samples  are 
intended  as  demonstrations  only.  Because  the  vessel  has  to  be  slit  and  laid  flat 
the  character  of  the  tissue  is  basically  altered.  The  tissue  samples  shown  in 
Figures  3.9  and  3.10  demonstrate  how  difficult  it  was  to  maintain  a  reasonably 
flat  surface  to  scan.  Figure  3.10  shows  tissue  that  is  relatively  normal  with 
only  a  small  portion  near  the  edge  with  some  calcification.  Figure  3.11  shows 
an  extreme  case  of  calcified  lesion  over  the  entire  surface. 


s-r/m  Assocf?tf<?rt  1 

CAST  SURFACE  COMPARATOR  i 


Figure  3.9.  Cast  surface  comparator. 


Figure  3.10.  Nearly  normal  tissue. 


44 


The  third  sample  does  not  have  a  photograph.  The  sample  had  a  rather 
unusual  structure  of  two  tubules,  perhaps  blood  vessels  or  scar  tissue  that 
formed  around  some  surgical  device,  that  crossed  the  main  vessel  section. 
This  portion  was  removed  when  the  vessel  was  slit  for  scanning.  The  rest  of 
the  vessel  is  relatively  normal. 


Reflection  Coefficient 


With  the  use  of  any  substitute  for  the  real  material  some  factors  can 
not  be  duplicated.  This  section  covers  the  differences  in  material  constants 
that  affect  the  reflection  of  acoustic  energy  from  a  surface.  The  reflection 
and  transmission  of  waves  from  an  interface  sure  determined  by  the  material 
characteristics  of  the  two  media. 

The  theoretical  basis  for  reflection  would  classify  the  reflection  from  the 
vessel  wall  to  be  more  like  a  liquid-liquid  interface  than  a  liquid-solid  interface. 
The  major  constituent  of  most  biological  tissue  is  water.  This  is  the  case  for 
both  blood  and  blood  vessels.  The  sample  surfaces,  however,  are  all  solids. 
The  reflection  from  them  will  vary  depending  on  their  physical  attributes  of 
p ,  density,  c,  sound  speed;  both  longitudinal  (cl)  and  shear  (cr),  and  these 
in  turn  depend  on  /?,  the  bulk  modulus  and  G,the  shear  modulus.  Table  3.4 
lists  the  relevant  values  (Kinsler  et  ad.,  1982,  pp.  461-462)(Shung,  1985,  p. 
310). 

The  reflection  and  transmission  coefficients  for  liquid-liquid  interfaces  are 
developed  in  the  following  manner  (Brekhovskikh,  1960,  pp.  16-18)  (Kinsler 
et  ad.,  1982,  pp.  131-133).  The  geometry  is  illustrated  in  Figure  3.12. 
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Table  3.4.  Material  constants. 


Material 

Sound 

Speed 

m/s 

Density 

kg/m 3 

Impedance 

(io-«) 

kg/mrs 

Bulk 

Modulus 

(*10“10) 

Pa 

Shear 
Modulus 
(•1(T10) 
Pa  ■  s/m 

■t!M 

“TO 

1.48 

M£iEM 

- 

mnom 

■MM 

1.65 

- 

- 

mniAffimw 

KM 

^■jqMNB 

7.6 

2.4 

m  mm 

mbim 

Kil 

ITS 

3.9 

2.5 

Lucite 

mmum 

'HI 

3.2 

MiBf.Mi 

Z 


Figure  3.12.  Reflection  geometry. 


The  incident,  reflected  and  transmitted  pressures  at  the  interface  are  de¬ 
fined  as: 

p.  _  p.giM-i^cosflj+tjzsinflj) 

p  _  p  e«(wi+ti*cos0r+fcixs»n#r) 


pt  _  -p  te*(u>t-k2*  cot  0t+k2xtmBt) 


(3.19) 


The  transmission  angle  is  assumed  complex  initially.  The  boundary  condition 
of  pressure  continuity  at  the  interface  requires  that 


_|_  p^ei(w<+*ix«n0r)  _  p<ci(u>t-fc2*«a*t) 


(3.20) 


For  this  to  be  true  for  all  x,  0,  =  6r-  It  also  follows  that  with  constant 
frequency, 


sin  <9, 
ci 


sin  ft 
C2 


(3.21) 


which  is  the  familiar  Snell’s  law.  The  magnitudes  of  the  functions  mean  that 


1  +  R  =  T, 


(3.2 2; 


where  the  reflection  ratio  is  defined  as  R  =  Pr/P,  and  the  transmission 
ratio  is  defined  as  R  =  P</P,.  The  second  boundary  condition  requires 
that  particle  velocity  normal  to  the  surface  is  zero.  This  results  in  a  second 
equation  in  R  and  T. 

1  -  R  =  (3.23) 

Z2  cos  ft 

where  z,  =  pm,  acoustic  impedance.  From  equations  3.22  and  3.23, 


R  _  (zj/zi)  ~  (cos  ft/  cos  9j) 
(z2/zi)  +  (cos  0t/  cos  ft ) 


(3.24) 
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The  value  of  Ot  can  become  complex  if  c\  <  ci  as  can  be  seen  from 

cos  $t  =  \f\  —  sin2  Ot  =  \/l  —  (C2/C1)2  sin2  0,  (3.25) 

when  (C2/C1)2  sin2  0,  >  1.  The  incident  angle  at  which  this  occurs  is  called 
the  critical  angle  and  is  defined  as 

sin  0C  =  —  (3.26) 

While  reflected  waves  are  of  primary  concern  for  this  thesis,  the  transmitted 
wave  information  is  included  for  completeness.  The  transmitted  pressure  then 
becomes 

pt  =  P  te-Te*ut~k*x  “  (3.27) 


where 

7  =  *2\/(c2/ci)2  sin20i  -  1 

The  transmitted  wave  travels  along  the  x  axis  and  decays  away  from  the 
boundary.  For  6t  real: 


_ 2*2/*! _ 

(*2/21)  +  (cos  6t/ cos  9i) 


(3.28) 


For  Ot  imaginary,  T  =  0.  The  reflection  power  level,  (201og(R)),  as  a  function 
of  incidence  angle  is  plotted  in  Figure  3.13. 
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Incident  Angle  (Degrees) 


Figure  3.13.  Reflection  power  level  for  aorta  tissue. 


The  theoretical  basis  for  the  liquid-solid  interface  is  similar  in  develop¬ 
ment.  The  result  is  somewhat  more  complicated  by  the  addition  of  shear 
waves  in  the  solid  (Brekhovskikh,  1960,  pp.  28-35).  The  particle  velocity  is 
defined  as 

v  =  grad<j>  -I-  curlxf;.  (3.29) 

where  <f>  is  the  longitudinal  potential  and  rp  is  the  transverse  (shear)  potential. 
The  potential  t/>  is  chosen  so  that  only  its  y  component  is  non-zero,  since  only 
the  x  and  z  components  will  contribute  due  to  the  orientation  selected.  The 
velocity  components  are  as  follows. 
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The  associated  wave  equations  are: 


V2<j> 


i  d2<t> 
4  dt2' 


V20 


1  d20 

4  dt 2 


(3.31) 


where  C£  and  cj  are  the  longitudinal  and  shear  velocities.  The  velocities  can 
be  expressed  in  terms  of  Lame’s  parameters  £  and  ft  or  equivalently  in  terms 
of  the  bulk  modulus,  /?  and  the  shear  modulus,  G. 


_  HHt 

CL  V  p 


(3.32) 


CL  = 


(3.33) 


Note  that  the  shear  velocity  is  always  less  than  the  longitudinal  velocity.  The 
stress  tensor  components  for  this  geometry  are: 


“  M  dx  +  dz>  +  *  dz’ 


r?  ,dux  dux 

z‘  =  "<  aT  +  a*  > 


(3.34) 


For  the  liquid  medium  t/>  =  0  and  n  =  0.  The  boundary  conditions  are 
continuity  of  Zz\ 


Aj  V2^i  =  A2  V2^2  4-  2W( 


92<h 

dz2 


+ 


dxdzh 


Zx  equal  to  zero: 

nd2<f>\  d2rl>i  _  d2ip2  _  n. 

dxdz  dx 2  dz 2 


d<j>\  _  d<h  drh 

dz  dz  dx 


and  continuity  of  uz: 


(3.35) 
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The  sound  wave  is  defined  by  the  following  formulas  for  the  incident  wave, 
<t>i ,  the  reflected  wave,  4>r,  the  transmitted  longitudinal  wave,  <fo,  and  the 


transmitted  transverse  wave,  02- 

<t>i  =  Ae*ui~k  (3.36) 

<t>r  =  ARe*(“t+k  (3.37) 

<t>2  =  ATe'(ut-k2'  co*  ““  (3.38) 

02  =  ASe,(w<-“2*co*7t+"2I,“  7t>  (3.39) 


The  wave  numbers  are  Arj,  and  /e2.  From  the  boundary  conditions  the  fol¬ 
lowing  relationships  are  found  for  the  reflection  and  transmission  coefficients 
using 

A:i  sin  0,  =  k 2  sin  =  k2  sin  7*  (3.40) 


and 


yields 


Zi  = 


i2  = 


T  = 


5  = 


Pi  ci 
cos  0, 


Z2  = 


P2C12 
cos  0t 


Zt  = 


P2CT2 

COS71 


Z2  cos2(27<)  +  Z<  sin2(27<)  —  Z\ 
Z2  cos2(27 1)  +  Z*  sin2(27<)  +  Zi 

_ 2Z2  cos(27t) _ 

Z2  cos2(27<)  +  Zt  sin2(27«)  +  Z\ 

2  Zt  sin  (271) _ 


(3.41) 

(3.42) 
(3.42) 
(3.42) 


Z2  cos2(27<)  +  Zt  sin2(27<)  +  Z\ 

Given  these  coefficients  there  are  three  regimes  of  operation.  When  0  < 
sin0,  <  ci/c£2>  both  6t  and  7 1  will  be  real  angles  and  the  reflection  coefficient 
will  be  real.  When  ci/c£2  <  sin  6,  <  ci/cn,  7<  will  be  real,  9t  will  be  complex, 
and  the  reflection  coefficient  will  be  complex.  The  longitudinal  transmitted 
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wave  will  travel  along  the  boundary  in  the  solid  and  the  shear  wave  will  be  a 
plane  wave.  When  c\/ct2  <  sin  0,,  both  $t  and  7 1  will  be  complex  angles  and 
the  reflection  coefficient  will  be  unity.  Both  transmitted  waves  will  propagate 
along  the  boundary.  Table  3.5  documents  the  critical  angles,  intensity  loss  at 
0  degrees  and  15  degrees  for  each  of  the  materials  and  aorta  tissue.  For  the 
three  materials  in  use  the  plots  for  their  reflection  power  level,  (201og(R)), 
are  shown  in  Figures  3.14,  3.15,  and  3.16  for  aluminum,  glass,  and  lucite, 
respectively. 


Table  3.5.  Critical  angles  and  intensity  loss. 


Material 

Critical 

An|les 

Intensity 

LossjdB) 

Intensity 

Loss(dB) 

15° 

Aorta 

- 

-21.6 

H0MM  III  1  ■  1 1 1 U 1 11 

-1-5 

Glass 

. -rs 

Lucite 

34.1 

-8.8 

^5 

5  10  15  20  25  30  35  40  45  50 


Incident  Angle  (Degrees) 

Figure  3.14.  Reflection  power  level  for  aluminum. 


5  10  15  20  25  30  35  40  45 


Incident  Angle  (Degrees) 
Figure  3.15.  Reflection  power  level  for  glass. 
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Chapter  4 


METHODS 


The  experiments  for  this  research  were  conducted  in  the  Ultrasound  Lab¬ 
oratory  of  the  Bioengineering  Department.  The  experiments  consist  of  two 
types.  The  first  experiment  type  is  used  to  plot  scattering  level  from  various 
surfaces  as  the  angle  of  incidence  is  varied.  This  is  used  to  establish  a  base¬ 
line  angle  for  investigating  the  backscatter  out  of  the  main  specular  lobe.  The 
second  type  of  experiment  examines  a  surface  at  two  angles.  The  first  angle 
is  at  0  degrees  from  the  vertical  to  the  mean  plane  of  the  surface.  The  second 
angle  is  at  15  degrees  from  the  vertical.  This  chapter  discusses  the  equipment 
used  and  the  procedure  followed. 


Experiment  Type  1 


The  purpose  of  the  first  experiment  was  to  find  the  amount  of  scattered 
energy  returned  from  a  smooth  surface,  various  sanded  surfaces  and  one  glass 
bead  surface  over  a  range  of  incident  angles.  The  system  has  a  monostatic 
transducer.  This  means  that  transducer  used  for  the  transmission  of  the 
pulse  is  also  used  as  the  receiver.  As  the  sample  is  rotated  the  incident  and 
reflection  angles  are  equal. 

The  procedure  for  this  experiment  is  to  mount  the  plexiglas  block  on  its 
adaptor  and  place  into  the  shaft  of  the  angular  measurement  vernier.  The 
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sanded  plexiglas  surfaces  are  oriented  so  that  their  parallel  sanded  directions 
are  perpendicular  to  the  incident  beam.  This  is  insures  that  the  beam  will 
see  the  roughest  surface  as  the  sample  is  rotated.  The  shaft  is  the  nearest 
vertical  shaft  in  the  figure.  The  vernier  is  the  squat  round  object  the  shaft  runs 
through.  The  sample  is  submerged  below  the  vernier.  The  second  vertical 
shaft  has  the  transducer  mounted  on  it.  The  beam  from  the  transducer  is 
aligned  to  get  the  maximum  return  from  a  stationary  object  on  the  vernier 
shaft.  After  this  alignment  the  pulse  is  bounced  off  the  surface  of  a  sample. 
The  sample  is  also  adjusted  to  get  the  maximum  return  level  at  specular 
reflection.  Once  a  peak  is  found  the  angle  is  recorded  for  use  as  the  starting 
point. 

The  first  surface  is  the  smooth  plexiglas  sample  for  use  as  a  baseline 
standard.  The  signal  return  is  routed  to  the  oscilloscope  and  the  peak  level  is 
aligned  with  a  point  on  the  screen.  Because  the  level  at  specular  reflection  is 
very  high  it  is  attenuated  through  a  series  of  calibrated  adjustable  attenuators. 
As  the  angle  is  changed  to  a  new  value  the  attenuators  are  set  to  adjust  the 
received  level  back  to  the  established  peak  value  on  the  oscilloscope. 

In  this  manner  each  surface  is  tested  and  the  results  recorded.  For  each 
the  starting  setting  is  the  same  as  for  the  smooth  surface.  In  this  manner  all 
surfaces  can  be  compared.  For  the  glass  bead  surface,  because  the  return  is 
very  diffuse,  the  smooth  reverse  side  is  set  up  and  then  the  vernier  is  rotated 
180°  to  arrive  at  0°  for  the  beads. 

For  this  first  experiment  set  the  following  equipment  is  used.  The  descrip¬ 
tion  of  each  piece  is  included  for  completeness  and  for  questions  of  repeatabil¬ 
ity.  The  transducer  is  driven  by  a  7.5  Megahertz  sine  wave  signal  generated 
with  a  Tektronix  FG504  Function  Generator.  The  signal  is  converted  to  a 
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pulse  with  a  Tektronix  PG505  Pulse  Generator.  The  pulse  is  5  microseconds 
long  and  has  a  pulse  repetition  rate  of  1  kilohertz.  The  rise  and  fall  rates  of 
the  pulse  are  1  microsecond.  The  attenuators  are  Kay  Elemetrics  Corporation 
Models  1/432D  and  432D.  The  Oscilloscope  is  a  Tektronix  SC  502.  The  flow 
chart  of  the  apparatus  for  this  experiment  is  shown  in  Figure  4.1. 


Figure  4.1.  Equipment  flow  chart  for  first  experiment. 


Experiment  Type  2 


The  second  experimental  type  is  designed  to  make  complete  C-mode  scans 
of  the  sample  surfaces  at  two  angles,  0°  and  15°  .  A  C-mode  scan  has  the 
azimuth  position  displayed  on  the  horizontal  axis  and  the  elevation  position 
displayed  on  the  vertical  axis  (IEEE  Standard  Radar  Definitions,  1982,  p.  9). 
The  view  is  the  same  as  viewing  the  sample  with  one’s  eye.  The  angles  chosen 
are  to  demonstrate  that  detail  lost  in  the  specular  (0°  )  scan  can  be  recovered 
in  an  oblique  (15°  )  scan.  Because  the  scan  rates  of  the  device  are  very  slow 
the  data  is  stored  in  a  computer  file  for  both  later  viewing  and  data  analysis. 
The  majority  of  the  scans  use  a  PDP-11  computer  with  a  Peritek  monitor 
for  displaying  the  scans.  The  scan  made  of  the  aluminum  cast  comparator  is 
done  with  a  Compaq  386. 

For  the  second  series  of  experiments  the  following  procedures  were  used. 
The  samples  were  set  in  position  in  the  same  manner  as  before  for  the  0°  scan. 
The  scan  is  controlled  by  the  computer.  The  computer  now  controls  the 
position  of  the  transducer,  by  the  use  of  three  stepper  motors.  The  stepper 
motors  can  be  seen  in  Figure  4.1  as  the  three  cylindrical  objects  oriented  in  the 
three  axes  directions  of  x  (transverse),  y  (vertical),  and  z  (longitudinal)  from 
the  alignm  'nt  for  the  transducer.  The  transducer  is  set  in  position  manually 
for  its  starting  position.  The  program  starts  and  the  computer  then  moves 
the  transducer  in  a  linear  raster  pattern  to  completely  scan  a  region  of  the 
sample.  The  width  and  height  of  the  scan  are  variables  for  the  program. 

The  level  is  stored  in  a  file  and  also  displayed  on  the  screen.  The  step  size 
for  the  PDP-11  scans  was  set  to  100pm.  The  Peritek  monitor  has  a  512  x  512 
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pixel  display.  This  can  display  a  region  as  large  as  the  samples.  The  scans  of 
this  size  unfortunately  can  take  up  to  two  days  to  complete.  The  PC  scans  are 
also  designed  to  be  512  points  wide.  Because  the  aluminum  cast  comparator 
is  11  centimeters  wide  compared  to  the  5  centimeter  width  of  the  blocks,  a 
larger  step  size  of  204 ftm  is  required  for  coverage.  Before  a  scam  is  conducted 
a  preview  scan  is  conducted  to  find  the  highest  likely  level  to  be  encountered. 
This  adlows  an  attenuation  level  to  be  set  to  keep  the  dynamic  range  of  the 
surface  within  the  stored  level  range  of  0  to  255. 

Once  the  scan  for  0°  is  made  the  transducer  returns  to  its  starting  posi¬ 
tion.  The  sample  is  rotated  to  15°  and  the  transducer  is  adjusted  in  x  and 
z  to  compensate  for  the  change  in  relative  position  to  the  sample.  A  second 
preview  scan  is  made  to  find  the  range  of  vadues  before  a  full  second  scan  is 
done. 

For  the  second  series  of  experiments  the  following  equipment  wais  used. 
The  transducer  is  a  Sound  Technologies  transducer  with  a  resonance  frequency 
of  10  megahertz.  The  equipment  is  adso  shown  in  Figure  4.2.  A  KB  Aerotech 
UTA-4  combines  many  of  the  functions  of  severed  devices  into  one.  This 
device  now  controls  the  pulse  level,  pulse  repetition  rate,  gated  return  signed 
and  attenuation  level.  A  Sony/Tektronix  390AD  receives  the  return  signed 
from  the  KB  Aerotech  and  digitizes  it  into  1024  points  of  ±  voltage  selected 
on  the  390AD  at  a  rate  of  60  megahertz.  The  digitized  values  are  then  sent 
to  the  computer  which  takes  the  absolute  value  of  the  voltage  eind  shifts  the 
value  right  by  one  bit  to  arrive  at  the  data  range  of  0  to  255.  The  flow  chart 
of  the  apparatus  for  this  experiment  is  shown  in  Figure  4.2. 


Figure  4.2.  Equipment  flow  chart  for  second  experiment. 


Chapter  5 


RESULTS 


This  chapter  presents  the  results  of  the  experiments  discussed  in  Chapter 
4  with  the  application  of  the  theory  presented  in  Chapter  2.  The  additional 
factor  of  differences  in  reflection  coefficients  presented  in  Chapter  3  is  also 
considered.  The  first  results  to  discuss  are  for  the  first  type  of  experiment. 
Next  the  results  for  the  second  experiment  set  are  shown.  The  second  set 
results  are  presented  in  two  forms.  The  first  is  the  qualitative  results  of 
the  plots  or  C-mode  images  of  the  scans.  The  second  form  is  the  statistical 
nature  of  the  levels  and  how  that  relates  to  the  theoretical  levels  expected. 
Finally  there  is  a  presentation  of  the  conclusions  of  this  thesis,  along  with 
some  suggestions  for  future  research. 


First  Experiment  Results 


The  first  set  of  experiments  show  the  level  of  the  backscattered  energy  for 
a  smooth  surface, the  four  sanded  surfaces,  and  a  beaded  surface.  The  levels 
are  tabulated  in  Table  5.1  and  are  referenced  to  the  reflected  level  at  normal 
incidence  for  the  smooth  surface.  The  angles  are  rounded  to  the  nearest  whole 
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Table  5.1.  (continued) 


Angle 

O 

Smooth 

Beads 

60 

Grit 

lOO 

Grit 

ISO 

Grit 

280 

Grit 

300 

Grit 

“31 

63.6 

2L2 

25.0 

“352 

63.6 

“5315 

“"5315 

42 

63.6 

20.6 

24.0 

37.6 

63.6 

63.6 

63.6 

43 

63.6 

20.6 

25.0 

36.4 

63.6 

63.6 

63.6 

44 

63.6 

21.3 

25.0 

36.6 

63.6 

63.6 

63.6 

45 

63.6 

22.0 

25.0 

37.4 

63.6 

63.6 

63.6 

The  value  for  I0  in  the  calculated  examples  uses  the  actual  values  from  the 
measurements  for  the  smooth  surface  sample.  Figure  5.1  shows  the  smooth 
surface  backscatter  referenced  to  the  specular  maximum. 

Now  that  an  established  I0  exists  the  equations  2.49,  2.50,  2.51,  and 
2.52  can  be  evaluated  and  compared  to  the  observed  levels.  Each  surface  is 
compared  in  the  following  plots.  Each  pair  of  plots  show  the  observed  levels 
in  the  first  plot  and  the  predicted  levels  are  shown  in  the  next  plot.  Figure  5.2 
shows  the  scattered  levels  from  the  400  grit  sanded  surface.  Figure  5.3  shows 
the  predicted  scattered  levels  from  the  400  grit  using  equations  2.49  and  2.50 
which  are  appropriate  for  g  =  0.0024.  Figures  5.4  through  5.11  show  results 
for  280,  180,  100,  and  60  grit  sanded  surfaces  and  the  predicted  levels  with  g 
values  of  0.0051,  0.0101,  0.0281,  and  0.0860,  respectively. 

Figure  5.12  shows  the  observed  levels  from  the  237 fim  diameter  bead  sur¬ 
face.  Figure  5.13  shows  the  predicted  scattering  level  range  from  the  237 /xm 
diameter  bead  surface  using  equations  2.49  and  2.51  for  g  =  1.17. 

The  plots  for  the  sanded  surfaces  of  100  and  60  grit  do  not  match  the 
level  of  the  predicted  values,  while  the  shape  does  show  good  agreement  with 
the  increased  level  in  the  vicinity  of  the  critical  angle.  This  would  imply  that 
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Figure  5.1.  Smooth  surface  backscatter. 


one  of  the  variables  in  equation  2.50  is  incorrect.  The  only  one  that  has  some 
question  about  it  is  the  correlation  length. 

The  correlation  length  required  to  match  the  level  for  the  100  grit  sample 
to  compensate  for  the  20  dB  higher  observed  level  would  be  95 pm,  which  is 
10  times  the  quantity  found  for  the  correlation  length.  The  correlation  length 
required  to  match  the  level  for  the  60  grit  sample  to  compensate  for  the  25 
dB  higher  observed  level  would  be  226/xm,  which  is  17.8  times  the  quantity 
extrapolated  for  the  correlation  length.  These  values  are  approximately  55 


Level  (dB) 


Figure  5.7.  Predicted  leveb  from  180  sanded  surface. 


Level  (dB) 


times  the  size  of  R a-  Using  similar  increases  for  the  other  sanded  surfaces 
causes  higher  than  observed  levels  to  be  predicted. 

For  the  case  of  the  beaded  surface  the  correlation  length  is  derived  from 
several  assumptions.  The  first  of  which  is  a  single  layer  of  beads.  For  this 
sample  the  assumption  is  probably  false  because  of  a  problem  in  making  the 
surface.  Once  this  assumption  is  dropped  the  value  of  both  Ra  and  X0  become 
unknown.  The  observed  levels  for  this  beaded  surface  are  10  dB  above  the 
higher  range  expected.  This  requires  an  increase  of  3.16  times  that  ratio  for 
a  similar  result. 

Second  Experiment  Results 

The  second  set  of  experimental  scans  has  the  advantage  of  providing  a 
visual  component  as  well  as  intensity  levels  to  evaluate  how  well  ultrasound 
can  characterize  roughness.  The  scans  come  from  two  computer  systems. 
The  first  system  is  a  PDP-11  with  a  Peritek  display.  The  Peritek  display 
has  512  x  512  pixels  and  can  display  256  levels  of  gray  scale  or  color.  All  of 
the  scans  except  the  aluminum  cast  comparator  are  made  with  the  PDP-11. 
This  computer  however  has  become  inoperable.  The  final  scans  are  done  on 
a  Compaq  386.  The  loss  of  the  PDP-11  requires  a  program  that  can  convert 
the  PDP  files  into  IBM  compatible  PC  files.  This  program,  which  allows 
the  manipulation  of  the  data  for  statistical  analysis,  is  documented  in  the 
Appendix. 

Over  fifty  scans  have  been  made  and  of  those  forty  are  available  for  analy¬ 
sis.  Of  these  eighteen  will  be  presented.  Table  5.2  lists  the  scans  by  file  name. 
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The  file  names  all  begin  with  the  letter  ‘D’  and  a  sequence  of  numbers  that  is 
the  date  of  the  scan  combined  with  the  incident  angle  of  the  scan.  The  only 
exceptions  are  the  scans  with  a  ‘C’  or  ‘N’  after  the  ‘D’  which  represent  the 
presence  of  carbon  particles  or  no  carbon  particles.  The  scans  are  listed  in 
Table  5.2. 


Table  5.2.  Scan  data  files. 


File 

Name 

WBSSSM 

Peak 

Level 

dB 

Description 

■imitM 

KWEK »!'■ 

wwm 

24 

Beaded  Strips 

IlM'll 

wwm 

38 

Beaded  Strips 

■M 

46 

Resolution  Pattern 

MM 

wnnm 

Resolution  Pattern 

■cHHEKiM 

mzam 

'  ~~  168231 

31 

■■MM 

■fi'OCUi 

mam 

■»miM 

iMK'l 

mm 

Kmi 

£87055 I 

■HiWililf 

imia 

wwm 

« ,tyyi  R* 

mmum 

wwm 

Km 

mamm 

m'Efi'l'M 

m>zvm 

Km 

088 

KMlKWiIifli 

T4 

“  088 

■  »I1W 

mmamm 

M>*zm 

kisih 

Carbon  Particles 

IflflBiHI 

wwm 

IKJ^I 

Carbon  Particles 

■»KIW 

m*um 

KiflH 

No  Carbon  Particles 

muaai 

timi 

wzvm 

mm 

No  Carbon  Particles 

K^KI 

Aluminum  Cast  Comparator 

Km 

Aluminum  Cast  Comparator 

The  following  pages  present  the  visual  evidence  in  a  16-step  decibel  gray 
scale.  The  gray  scale  is  illustrated  in  Figure  5.14.  The  lowest  level,  0,  is  at 
the  left  of  the  figure.  The  levels  are  as  represented  in  Table  5.3.  The  dynamic 
range  of  a  256  level  system  measuring  voltage  is  48.1  dB.  The  colormap  has 
a  maximum  of  64  levels  in  a  range  0  —  63.  The  log  display  scale  shifts  the 
emphasis  on  higher  levels  and  the  variation  between  those  levels. 


Level 

Point 

Palette 

Peak 

Range 

Level 

dB 
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The  scans  are  presented  as  pairs,  with  the  0°  scan  first,  followed  by  the 
15°  scan.  These  are  all  shown  at  approximately  twice  the  true  scale  of  the 
sample.  The  exception  is  the  aluminum  cast  comparator  scans,  which  are 
1.25  :  1  scale.  The  scans  are  presented  in  roughly  chronological  order  and  will 
be  discussed  in  this  order. 

The  first  set  of  scans  is  on  the  resolution  block  shown  in  Figures  5.15  and 
5.16.  The  images  includes  the  central  10  millimeter  square  and  some  portions 
of  each  of  the  other  smaller  squares.  This  result  is  encouraging,  demonstrating 
that  surface  roughness  differences  '•an  be  observed.  This  also  demonstrates 
that  for  extremes  in  roughness  even  the  normal  scan  can  show  differences  in 
surface  roughness. 

The  next  two  sets  of  scans  are  of  the  same  surface,  the  beaded  strips.  The 
first  set,  D52115  and  D52200,  are  of  the  same  area  however  the  0°  scan  is  only 
half  the  size  because  of  the  time  required  to  make  the  larger  scan  and  poor 
scans  that  resulted  due  to  computer,  stepper  motor,  and  transducer  problems 
with  long  scan  times  (longer  than  24  hours).  The  larger  scan  does  overlap  the 
region  covered  in  the  small  scan.  Several  features  can  be  seen.  The  outer  two 
rows  of  larger  beads  axe  not  in  these  seems.  In  the  0  degree  scan  the  larger  two 
sizes  of  beads  in  the  center  scatter  more  of  the  energy  away  from  the  specular 
direction  than  the  three  rows  of  smaller  beads  on  either  side.  In  the  15°  scan 
only  one  of  the  larger  bead  rows  appears  to  make  a  larger  contribution  to  the 
reflected  angle  energy. 


Figure  5.15.  Resolution  block  0  degrees. 


Figure  5.16.  Resolution  block  15  degrees. 


Figure  5.17.  D52200;  bead  strips  0  degrees. 


Figure  5.18.  D52115;  bead  strips  15  degrees. 
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The  second  set  of  bead  strips  does  include  portions  of  the  ten  strips.  The 
scans  are  shown  in  Figures  5.19  and  5.20.  Again  the  number  of  scan  lines  was 
reduced  for  time  considerations.  In  Figure  5.19  the  larger  beads  in  the  center 
and  at  the  edges  have  a  lower  return  than  the  other  beads  in  the  0°  scan.  The 
reverse  is  true  for  the  15°  scan. 


Figure  5.19.  D52900;  bead  strips  0  degrees. 


Figure  5.20.  D52915;  bead  strips  15  degrees. 


The  next  two  sets  of  scans  are  of  plexiglas  blocks  with  regions  that  have 
glue  applied.  In  the  first  surface,  two  areas  of  carbon  particles  sprinkled  on 
glue  are  separated  by  a  sanded  region  (280  grit).  The  sanded  region  is  a 
central  strip  in  the  middle  of  the  block  from  top  to  bottom  and  across  the  top 
of  the  right  hand  side.  The  area  to  the  left  side  has  particles  of  10/xm.  The 
area  to  the  right  has  particles  of  60 pm. 

The  second  surface  has  only  the  glue  on  the  two  areas  separated  by  a 
sanded  region.  The  carbon  particle  scans  are  shown  in  Figures  5.21  and  5.22. 
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The  glue  only  scans  are  shown  in  Figures  5.23  and  5.24.  There  are  some 
interesting  regions  that  can  be  seen  in  these  figures.  In  both  of  the  0°  scans 
the  region  of  sanded  surface  appears  uniform.  The  particles  on  glue  region 
and  the  glue  only  region  have  a  similar  undulating  appearance.  The  15°  scans 
also  look  quite  similar.  Unfortunately  the  regions  of  highest  return  do  not 
seem  to  be  aligned  with  either  of  the  regions  of  particles  on  glue  or  glue  only. 
Another  feature  that  can  be  seen  is  many  small  dark  dots.  These  are  pixels 
with  very  low  levels  (<  12).  They  appear  to  be  data  drop  out  problems  with 
the  digitizer  board  on  the  PDP-11.  They  appear  uniformly  throughout  all 
four  scans.  The  worst  case  is  in  Figure  5.22,  which  is  also  the  last  of  the  four 
to  be  made. 

The  next  set  of  scans  is  of  the  aluminum  cast  comparator.  The  scans  are 
shown  in  Figures  5.25  and  5.26.  The  roughest  surface  is  on  the  left.  The 
highest  levels  on  the  0°  scan  appear  to  be  from  the  central  region.  This  is  not 
the  expected  result.  Upon  further  investigation  the  plate  was  found  to  have 
a  slight  curvature  to  the  surface.  Because  of  the  sensitivity  of  the  scattering 
function  around  0°  the  slight  curvature  is  enough  to  cause  a  disruption  in 
the  expected  result.  Measurements  of  the  plate  were  made  to  determine  the 
curvature.  The  result  is  that  over  the  whole  plate  the  normal  to  the  surface 
vaxied  by  1°  to  2°  depending  on  where  the  measurements  are  taken.  The 
15°  scan  shows  somewhat  similar  problem  with  the  highest  levels  coming  from 
the  roughest  areas  (expected)  and  from  the  centred  region  again  (unexpected). 
Because  of  the  curvature  the  angle  of  incidence  is  varying  Mound  the  region 
of  the  critical  single  for  aluminum  of  13.6°.  Nonetheless  the  scan  does  show 
variations  for  each  surface  roughness  area. 
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Figure  5.21.  DC120;  carbon  particles  0  degrees. 


Figure  5.22.  DC1315;  carbon  particles  15  degrees. 


The  next  three  sets  of  scans  are  of  tissue  samples.  As  noted  before  in 
Chapter  3,  the  surfaces  retained  some  of  their  natural  curvature  when  open 
and  flat.  The  scans  here  are  demonstrations  of  the  feasibility.  They  show  that 
as  the  angle  is  changed  some  variation  will  occur  and  the  cause  of  the  changed 
level  is  the  variation  in  the  surface.  Unfortunately  most  of  the  variation  seen  in 
the  scans  is  due  to  the  undulation  from  a  curved  cylindrical  surface  reacting  to 
internal  stresses  rather  than  the  actual  roughness  of  the  surface  in  its  normal 
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Figure  5.23.  DN120;  no  carbon  particles  0  degrees. 


state.  Figures  5.27  and  5.28  show  a  relatively  normal  sample.  This  is  the  same 
sample  shown  in  Figure  3.10.  The  regions  of  highest  return  are  from  the  curve 
of  the  surface.  Figures  5.29  and  5.30  show  another  relatively  normal  sample. 
There  are  regions  here  that  have  a  large  number  of  data  dropouts  in  regions 
of  higher  return.  Again  the  regions  of  highest  return  are  from  the  curve  of 
the  surface.  Figures  5.31  and  5.32  are  the  scans  that  correspond  to  the  highly 
calcified  tissue  pictured  in  Figure  3.11.  Here  there  is  some  indication  by  more 


Figure  5.25.  D70700;  aluminum  cast  comparator  0  degrees. 


Figure  5.26.  D70915;  aluminum  cast  comparator  15  degrees. 

frequent  undulation  of  levels  that  the  plaque  is  being  shown  as  well  as  the 
tendency  of  the  surface  to  recurve  itself  back  into  a  cylinder. 

The  next  section  describes  some  statistical  relationships  found  by  using 
the  program  documented  in  the  Appendix.  The  regions  of  various  scans  se¬ 
lected  are  shown  again  with  the  areas  of  interest  drawn  as  rectangles  on 
the  scans.  These  levels  are  compared  with  the  plots  of  expected  scattering 
strength  as  in  the  first  experiment  set.  Because  the  source  of  disagreement  in 
the  first  experiment  set  appears  to  be  correlation  length,  an  attempt  is  made 
to  find  the  correlation  length  from  the  backscatter  field. 

For  each  of  the  scans  the  following  analysis  is  made.  A  section  of  a  scan  is 
selected  that  has  the  same  value  of  roughness.  The  data  points  are  averaged 
and  the  standard  deviation  found  for  that  area.  The  decibel  level  is  computed 
for  the  average  value,  the  average  value  plus  one  standard  deviation  and  the 
average  value  minus  one  standard  deviation.  The  decibel  level  is  referenced 


Figure  5.27.  D52700;  A87065  tissue  sample  0  degrees. 
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Figure  5.28.  D52715;  A87065  tissue  sample  15  degrees. 

to  the  absolute  peak  level  in  the  0°  scan.  The  assumption  is  that  the  highest 
level  is  equivalent  to  a  plane  surface  somewhere  in  the  scan.  Because  there 
is  a  difference  in  attenuation  for  the  15°  levels,  an  adjustment  is  made  by 
subtracting  the  difference, 

•^15  =  ^15, measured  (Ao  A15)  5.1 

where  Aq  is  the  attenuation  level  for  the  0°  scan  and  Ais  is  the  attenuation 
for  the  15°  scan.  In  addition,  for  each  rectangular  area  selected,  at  least 
five  lines  were  chosen  to  find  the  correlation  lengths  for  the  acoustic  field. 
Finally  the  three  averaged  levels  for  0°  and  15°  are  compared  with  a  plot  of 
the  the  predicted  backscatter  levels  using  the  calculated  or  derived  values  of 


Figure  5.29.  D71900;  088  tissue  sample  0  degrees. 


Figure  5.30.  D71915;  088  tissue  sample  15  degrees. 


the  correlation  length  set  in  the  first  experiment  set  and  with  predicted  levels 
using  the  average  correlation  length  from  the  acoustic  held. 

The  first  surface  scan  to  be  analyzed  is  the  resolution  square.  Figures  5.33 
and  5.34  show  the  scans  and  the  two  regions  used  to  compute  the  averages. 
The  regions  are  enclosed  in  rectangles.  The  average  level  and  the  standard 
deviation  of  the  levels  in  each  region  as  well  as  the  correlation  lengths  of  lines 
in  each  region  are  calculated  by  the  computer  program.  The  average  levels 
and  the  average  correlation  length  are  shown  in  Table  5.4.  The  roughness  of 


Figure  5.31.  D52600;  168231  tissue  sample  0  degrees. 


Figure  5.32.  D52615;  168231  tissue  sample  15  degrees. 


the  rough  area  is  unknown.  The  smooth  area  however  was  made  with  400 
grit  sandpaper.  Because  the  sanded  surface  has  a  value  g  much  less  than 
unity,  equations  2.49  and  2.50  are  used  for  the  prediction  curves.  Figure  5.34 
shows  how  the  levels  compare  with  the  correlation  values  from  Table  3.1.  In 
this  figure,  as  in  all  the  following  the  average  level  measured  is  marked  with 
a  circle  on  the  graph,  o.  If  more  than  one  area  is  measured  for  the  same 
surface  roughness,  there  is  a  circle  for  each  average.  In  addition  the  plus  and 
minus  one  standard  deviation  levels  are  marked  with  a  diamond,  o.  Thus  the 
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extent  of  the  markings  at  0  and  15  degrees  show  the  confidence  interval  for 
the  measured  values.  The  average  correlation  length  for  the  acoustic  field  was 
1340/im.  The  comparison  is  shown  in  Figure  5.35. 

Table  5.4.  Backscatter  levels  for  400  sanded  surface. 


The  other  scans  that  include  a  portion  of  sanded  surface  axe  the  four  with 
and  without  carbon  particles.  The  correlation  length  for  the  sanded  regions 
was  calculated  from  the  two  central  regions.  The  scan  areas  are  shown  in 
Figures  5.37  and  5.38.  Only  the  scans  with  the  carbon  particles  are  shown 
because  there  is  essentially  no  difference  between  the  carbon  and  carbonless 
scans.  Figure  5.39  shows  the  average  backscatter  levels  compared  with  the 
predicted  values  using  A0  from  the  Table  3.1.  The  average  levels  and  the 
standard  deviations  are  shown  in  Table  5.5.  The  sanded  area  was  divided 
into  three  segments  on  the  two  types  of  scans.  Therefore  there  are  a  total  of 
six  entries  for  each  angle.  The  values  are  presented  in  the  following  order.  The 
first  three  values  for  each  angle  are  from  the  scam  with  the  carbon  particles. 
The  next  three  are  from  the  scan  without  the  carbon  particles.  Of  the  three, 
the  first  is  from  the  upper  central  rectangle.  The  second  of  the  three  is 
from  the  lower  central  rectangle.  The  third  of  three  is  from  the  upper  right 
rectangle.  The  values  are  reasonably  close  to  one  another.  Figure  5.40  shows 
the  comparison  of  the  observed  backscatter  levels  compared  with  the  predicted 
level  with  the  acoustic  field  value  for  the  correlation  length  of  1385/im. 


Figure  5.33.  Resolution  block  0  degrees. 


Figure  5.34.  Resolution  block  15  degrees. 


The  results  for  the  beaded  strips  are  presented  in  Table  5.6  in  which  data 
are  presented  for  each  angle  a  maximum  of  four  times.  There  are  two  strips 
on  the  surface  and  because  there  are  two  sets  of  scans;  D52200  and  D52115 


Figure  5.37.  DC120;  carbon  particles  0  degrees. 


Figure  5.38.  DC1315;  carbon  particles  15  degrees. 


as  well  as  D52900  and  D52915.  The  values  vary  greatly  at  each  angle.  The 
first  set  has  only  eight  columns  as  can  be  seen  in  Figures  5.41  and  5.42.  The 
second  set  has  ten  areas  as  shown  in  Figures  5.43  and  5.44. 
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Figure  5.41.  D52200;  bead  strips  0  degrees. 


Figure  5.43.  D52900;  bead  strips  0  degrees. 


Figure  5.44.  D52915;  bead  strips  15  degrees. 


92 


The  predicted  levels  for  backseat  ter  level  from  the  beaded  surfaces  use 
equations  2.49  and  2.51  or  equation  2.52  alone,  depending  on  the  value  of  g. 
For  the  two  smallest  sizes  of  beads  equations  2.49  and  2.51  are  used.  For  the 
other  sizes  both  equations  2.49  and  2.51  as  well  as  2.52  are  used  and  compared. 
This  is  done  because  the  value  of  g  for  these  sizes  is  on  the  borderline  for  the 
conditions  governing  the  theory.  If  equations  2.49  and  2.51  are  used,  there  are 
two  lines  plotted  that  establish  the  predicted  range  of  backscatter  intensity. 
If  equation  2.52  is  used  there  is  only  a  single  line  shown.  Figures  5.45  and 
5.46  show  the  backscatter  levels  plotted  on  the  predicted  levels  for  the  50 fim 
beaded  surface  for  the  calculated  A0  from  Table  3.2  and  for  the  measured 
acoustic  X0,  respectively.  Figures  5.47  through  5.50  show  similar  plots  for 
the  105 fim  and  150/im  diameter  bead  surfaces  using  equations  2.49  and  2.51. 
Figures  5.51  and  5.52  show  plots  for  the  150^m  diameter  bead  surface  using 
equation  2.52  assuming  that  g  is  much  greater  than  1.  The  next  four  plots 
show  the  measured  data  for  the  215/im  diameter  bead  surface.  Figures  5.53 
and  5.54  show  the  observed  data  compared  to  equations  2.49  and  2.51.  Figures 
5.55  and  5.56  show  the  observed  data  again  using  equation  2.52.  Figures  5.57 
through  5.60  show  similar  plots  for  the  237 \im  diameter  bead  surface  using 
equations  2.49  and  2.51,  and  2.52. 

The  last  surface  to  be  shown  is  the  aluminum  cast  comparator.  The 
problem  with  the  curvature  has  already  been  discussed.  The  backscatter 
strength  levels  are  presented  in  Table  5.7.  The  spread  of  levels  across  the  one 
standard  deviation  is  much  tighter  than  for  the  beads.  The  trend  expected  of 
higher  levels  at  0°  incidence  angle  for  the  least  rough  is  not  seen.  The  reverse 
trend  for  the  15°  incidence  angle  is  also  not  consistently  present.  The  regions 
used  to  make  the  averages  are  shown  in  Figures  5.61  and  5.62.  Comparisons 
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Figure  5.47.  105 /xm  bead  surface  (g  ~  1)  with  material  A0. 
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Figure  5.48.  105 pm  bead  surface  (p  ~  1)  with  acoustic  A0. 
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Figure  5.56.  215fim  bead  surface  (<7  1)  with  acoustic  X0. 
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of  the  correlation  length  values  estimated  from  Table  3.3  and  the  measured 
values  using  the  acoustic  field  are  shown  in  Figures  5.63  through  5.84.  The 
first  three  surfaces  have  roughness  values  that  meet  the  following  condition: 
j  <  1.  These  are  illustrated  in  Figures  5.63  through  5.68.  The  next  four 
surfaces  meet  the  condition,  g  ~  1.  The  comparisons  for  these  surfaces  are 
shown  in  Figures  5.69  through  5.76.  The  two  curves  illustrate  the  upper  bound 
and  the  lower  bound.  For  the  last  two  surfaces  the  condition,  g  ~  1,  or  the 
condition,  g  <£.  1,  might  apply  so  both  are  shown.  These  are  in  Figures  5.75 
through  5.84.  For  some  of  the  acoustic  field  correlation  lengths  the  predicted 
field  has  a  gain  predicted  that  is  impossible.  This  casts  considerable  doubt  on 
the  relationship  between  the  measured  acoustic  field  correlation  length  and  its 
actual  value.  Table  5.8  compiles  a  list  of  all  the  roughness  of  the  surfaces  in 
this  study  along  with  the  correlation  lengths,  both  found  or  derived  and  the 
acoustic  field  correlation  length.  In  addition  the  table  gives  the  values  of  g  for 
all  the  materials.  Tables  5.9  and  5.10  compile  the  errors  for  all  the  surfaces  at 
0°  and  at  15°  ,  respectively.  The  levels  for  each  surface  type  take  the  average 
of  the  observed  levels  if  more  than  one  value  exists  for  a  particular  surface. 
For  the  conditions  where  equation  2.51  is  used,  the  error  is  only  listed  if  it  is 
outside  the  bounds  of  the  two  curves.  The  error  is  calculated  from  the  nearest 


line. 


Figure  5.61.  D70700;  aluminum  cast  comparator  0  degrees. 


Figure  5.62.  D70915;  aluminum  cast  comparator  15  degrees. 


Level  (dB)  Level  (dB) 


103 


Incident  Angle  (Degrees) 


Figure  5.64.  0.51  fimR^  aluminum  surface  3  <  1  with  acoustic  A„. 


104 


Incident  Angle  (Degrees) 


Figure  5.65.  1.53 nmR*  aluminum  surface  g  <  1  with  material  A0. 
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Figure  5.66.  1.53 pmR,A  aluminum  surface  g  <  1  with  acoustic  A„. 
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Figure  5.67.  3.05 iimR*  aluminum  surface  g  <1  with  material  X0. 


Figure  5.68.  3.05 nmR^  aluminum  surface  1  with  acoustic  X0. 
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Figure  5.69.  5.08 nmR^  aluminum  surface  g  ~  1  with  material  \0. 


Figure  5.70.  5.08 nmR^  aluminum  surface  g  ~  1  with  acoustic  A0. 
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Figure  5.72.  7.62 pimR^  aluminum  surface  g  ~  1  with  acoustic  A0. 
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5.73.  10.67 hthRa  aluminum  surface  g  ~  1  with  material  A0. 
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Figure  5.74.  10.67/imi?yi  aluminum  surface  g  ~  1  with  acoustic  A0. 


Level  (dB)  Level  (dB) 


109 


Incident  Angle  (Degrees) 


Figure  5.76.  14.22/xmi? &  aluminum  surface  g  ~  1  with  acoustic  X0. 
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Figure  5.77.  18.29 nmR^  aluminum  surface  g  ~  1  with  material  \0. 


Figure  5.78.  18.29 hthRa  aluminum  surface  g  ~  1  with  acoustic  A0. 
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rigure  5.79.  18.29{imR.A  aluminum  surface  .</  »  1  with  material  A„. 


Figure  5.80.  18.29/imfix  aluminum  surface  g  »  1  with  acoustic  X0. 
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5.81.  22. 86 pm  RA  aluminum  surface  g  ~  1  with  material  Ac. 


Figure  5.82.  22.86pmRA  aluminum  surface  g  ~  1  with  acoustic  X0. 
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Table  5.8.  Material  and  statistical  properties. 


Ra 

[Am 

9 

(0°) 

Ao 

(material) 

[im 

A* 

(acoustic) 

[Am 

Material 

0.51 

1.81  •  10~3 

2.8 

1220 

Aluminum 

0.58 

2J36  •  10“3 

4.1 

1340 

Lucite 

0.85 

5.07  •  10-'6 

5.0 

1385 

Lucite 

1.20 

10.11  •  10"3 

7.5 

- 

Lucite 

1.52 

16.22  - 10“3 

8.3 

1180 

Aluminum 

2.00 

28.07  10"3 

9.5 

- 

Lucite 

3.05 

65.29  10”3 

16.5 

1260 

Aluminum 

3.50 

85.98 -10"3 

12.7 

- 

Lucite 

5.08 

0.181 

27.43 

980 

Aluminum 

7.25 

1069 

T75T 

I960 

Glass 

7.62 

ras 

41.15 

TT80 

Aluminum 

10.67 

07799 

57762 

1400 

Aluminum 

14.22 

019 

76779  ~ 

T300 

Aluminum 

15.23 

1.628 

16.41 

1240 

Glass 

18.29 

2348 

95777 

1950  H 

Aluminum 

21.75 

3.320 

23745 

1270 

Glass 

22.86 

r  3.668 

123.44 

1300 

Aluminum 

31.18 

r  6323 

33.59 

1300 

Glass 

34.37 

57291 

1  37.03 

r  890  — 

Glass 
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Chapter  6 


CONCLUSIONS 


The  exercise  undertaken  shows  that  the  technique  has  the  ability  to  defi¬ 
nitely  differentiate  different  levels  of  roughness.  The  scan  images  are  the  best 
evidence  of  the  feasibility.  The  supporting  statistical  data  are  not  as  encour¬ 
aging.  This  is  not  too  surprising  considering  the  nature  of  Kirchhoff  theory. 
The  theory  is  only  exact  when  the  surface  is  a  smooth  plane  (Ogilvy,  1991, 
p.  100).  As  the  surface  becomes  more  and  more  rough,  errors  increase.  Still 
it  is  the  only  theory  that  provides  a  basis  for  comparison. 

Some,  if  not  all,  of  the  problems  seen  with  the  data  can  be  attributed 
to  two  problems.  The  first  is  the  uncertainty  of  the  correct  values  for  cor¬ 
relation  length.  The  various  attempts  of  measuring  the  value  have  not  been 
successful.  The  smallest  resolution  of  the  data  points  using  the  strip  chart 
yields  correlation  lengths  equal  to  the  sampling  interval.  The  attempt  to  use 
the  acoustic  field  has  improved  the  correlation  with  some  of  the  predicted 
backscatter  values.  However,  this  technique  also  has  some  clearly  false  results 
with  the  aluminum  cast  comparator. 

The  second  problem  is  the  model  surfaces.  The  manually  fabricated  sur¬ 
faces  using  the  plexiglas  blocks  prove  to  be  inconsistent.  The  primary  reason 
for  this  is  the  lack  of  any  standard  surfaces  in  the  range  of  surface  roughness 
of  interest.  The  one  manufactured  standard,  the  cast  surface  comparator, 
while  a  standard  for  roughness,  had  no  documentation  for  the  correlation 
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length  and  the  surface  curvature  causes  considerable  difficulty  for  measure¬ 
ment.  One  problem  with  the  bead  strips  also  can  not  be  ignored.  The  width 
of  the  strips  was  approximately  equal  to  the  beamwidth  of  the  transducer. 
This  was  not  a  judicious  choice  for  acquiring  a  large  number  of  points  from  a 
single  roughness  region. 

For  continued  research  in  this  area,  the  following  suggestions  are  appro¬ 
priate.  There  is  a  need  for  reliable  standard  surfaces  that  are  in  the  range 
of  surface  roughness  required.  There  is  also  the  need  for  research  into  the 
roughness  of  the  lesions  themselves.  Beyond  repeating  these  same  experi¬ 
ments,  the  next  step  would  be  the  use  of  a  rotating  head  transducer  similar 
to  the  commercial  intravascular  probe  on  the  interior  of  cylindrical  standard 
surfaces  and  in  vitro  tissue  samples. 
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Appendix 


SOFTWARE 


This  appendix  includes  a  listing  of  the  program  replot. cpp.  replot  is  a 
program  that  has  evolved  during  the  development  of  this  thesis  into  an  all 
inclusive  analysis  tool  for  the  scans  made  with  the  PDP-11  and  the  Compaq 
386.  When  the  scans  are  made  using  the  PDP-11  with  the  Peritek  display, 
the  data  are  stored  in  a  compressed  form  along  with  two  color  maps  and  a 
header.  To  continue  using  the  data  it  became  necessary  to  develop  a  program 
that  could  redisplay  the  already  collected  data  on  another  computer.  The 
data  files,  on  8”  diskettes,  were  taken  to  the  Electrical  Engineering  VAX  and 
transferred  to  5  diskettes  usable  in  personal  computers.  Considerable  effort 
was  expended  to  decode  the  compression  algorithm  used  by  the  Peritek  board 
to  again  display  the  data.  The  original  data  have  256  levels  of  color  or  gray 
scale.  The  VGA  display  most  commonly  used  on  personal  computers  has  16 
levels  of  color.  The  program  writes  a  file  of  uncompressed  data  and  then  shifts 
the  data  right  by  four  bits  and  then  displays  the  data.  The  original  Peritek 
display  is  a  512  x  512  display  while  a  VGA  display  is  640  wide  and  480  high. 
The  data  when  displayed  on  the  VGA  are  centered  horizontally  and  the  last 
line  of  data  is  aligned  with  the  bottom  of  the  display.  This  means  that  32 
lines  at  the  top  are  omitted.  This  was  chosen  because  the  original  scans  are 
at  the  bottom  if  they  sure  near  an  edge.  For  the  scans  made  later  the  data 
have  32  lines  of  blank  data  in  the  beginning  in  order  to  align  the  scan  of  the 
aluminum  cast  surface  comparator  with  the  top  of  the  display. 
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//  This  prograa  takas  a  fils  from  a  Paritak  program  and  // 

//  redisplays  it  11  tha  lila  already  exists  in  tha  torn  raquirad  // 

//  lor  display  on  a  PC  tha  prograa  aill  display  it  and  allow  // 

//  ragiona  o 1  tha  display  to  ba  salactad  with  a  aousa  and  tha  // 

//  aaan,  standard  deviation,  and  paak  walua  to  ba  lound.  It  also  // 

//  allows  lines,  both  vortical  and  horizontal  to  ba  salactad  and  // 

//  tha  corralatioa  langth  lound  along  thoaa  linos  // 

tinduda  < graphics . h>  //  For  graphics  library  lunctions  // 

tinduda  <stdlib.h>  //  For  axitO  // 

tincluda  <stdio.h> 

tinduda  <conio.h> 

tincluda  <alloc.h> 

tincluda  <string.h> 

tincluda  <aath.h> 

tincluda  <dos . h> 

tincluda  "aousa2.iH 

Hrasult*  Rasult; 

unsignad  int  plot .peak  -  0; 

Hstatus  Posit  ionJl,  Position!.; 

Host  db_ratio,  illoat,  awaraga,  average_2,  st.dev,  data.lloat,  n_ points. 
n.points.oll, cor .walua,  l.data,  data. array [1066] ,  corrarray [S12] ; 

//  Global  wariablss  —  sat  by  calc.coordsO  // 

int  aax.x,  aax.y ,  aid.x ,  aid.y ;  //  Haxiaua  x-  and  y-coordinatas  // 

unsignad  char  coap.dava.  rag_data,  aap.l ,  aap_2,  gllag,  data.byta; 

unsignad  char  color _aap_uppar_word  ■  0x0 F,  color_aap.lowar.word  ■  OxFF; 

unsignad  long  color .aap.l [256] .  tost .walua,  haad4,  haadS; 

unsignad  long  col or .map. 2 [256] ; 

unsignad  char  haadar.worda [46] ; 

unsignad  char  aap.lull[512] ,lina_scan[5l2] ; 

char  n_boxes[3]  ,xl[4]  ,  yl[4]  ,  x2[4]  ,  y2[4]  ; 

unsignad  int  dacoda .header  (unsigned  char  w[48] ,  unsignad  int  word.count) ; 

struct  palattatypa  lar  apal.pt,  spal.pt.oldalULL; 

struct  palattatypa  lar  pal,  pal  .old; 

void  stora.aap_l(void); 

void  stora.Bap_2(void); 

void  straight .writa (void); 

void  unpack (void); 

void  statistics(void) ; 

void  corral (void) ; 

void  linaar.gr ay  (struct  palattatypa  pal) ; 
void  log_gray(struct  palattatypa  pal) ; 
void  old.palatta(struct  palattatypa  pal); 
void  naw.palatta (struct  palattatypa  pal) ; 
void  joa.palatta(struct  palattatypa  pal) ; 


int  sat^raph(void) ;  //  Initializa  graphics  // 
void  calc.coords(void);  //  Scale  distances  onscreen  // 
void  draw.scan(void) ;  //  Draw  scan  // 
void  draw_rect(void);  //  Draw  rectangles  Iron  a  lila  // 
void  draw_zooa(void) ;  //  Draw  a  zoom  scan  // 


void  surlace.level(void) ;  //  Calculate  levels  in  a  region// 

void  corral.lina(void) ; 

void  gat .key (void) ;  //  Display  text  on  graphics  screen// 

//  wait  lor  key  // 

void  gat .ptl (void);  //  Find  lirst  point  // 

void  gat _pt2 (void) ;  //  Find  second  point  // 

int  peak.datact (unsigned  char  data_point);  //  Find  highest  level  in  plot  // 
unsignad  int  color.laval,  head [6]; 
unsigned  int  paak.to.now; 

int  loop  ct,  line.ct,  array.indl,  scan.x,  scan.y, 

EXIT  -  FALSE,  EXIT!  -  FALSE; 


int  start.x,  end.x,  start.y,  end.y,  i_a,  j.a,  index.cor; 
char  ana .  ansi,  b.char  *  ’ 

char  filename  [30]  ,  fil*name2[30] ,  outf  il«nasa[30]  ,  outf ilename2[30]  , 
filenames [30] ,  outf ilena»e3[30] ; 
char  inbuf  Cl] ; 

FILE  *1 ils.pt; 

FILE  *f ile.pt 1; 

FILE  *fils.pt2; 

FILE  *fils.pt3; 

FILE  sfils_pt4; 

int  uppsr_lsft_x,  upper.left.y,  loser .right.x,  loser.right.y ,  tssp.coor 

int  del.x.del.y,  npoint,  boxss,  n.box; 

int  dslta.x,  delta.y,  rsax.linit,  ysax.linit; 

int  TPos(  int  TP,  int  Los,  int  High  ) 

{ 

rstum(  (  TP  >-  Los  )  ft*  (  TP  <-  High  )  ); 

} 

to id  BoxItaaC  int  x,  int  y,  int  s,  int  h,  char*  text  ) 

{ 

ssttsxt justify (  CEITER.TEXT,  CEHTER.TEXT  ); 
rectangle (  x,y,x+s,y+h  ); 
outtsxtxy(  x*(e/2).  y-Kh/2)  .text  ); 

} 

int  sain  (sold) 

{ 

//  Find  out  shich  *.sch  fils  to  bs  plottsd  // 
clrscrC  ); 

print!  ("  Doss  ths  data  fils  already  exist  for  plotting  on  a  PC?  \n  " 

"  inssar  y  or  T  if  it  exists.  \n"); 
gets (inbuf) ; 

sscanf  (inbuf ,"Xlc",ftans); 
if  ((ana!-'y*)  U  (ana 

printf  ("Which  VCH  data  fils  do  you  sant  plotted?  ") ; 
scanf  ("Xs" ,f ilenaae) ; 

if  ((fils.pt  -  fopsn  (filename,  "rb"))  “  HULL) 

{ 

f printf  (stdsrr,  "Cannot  open  input  fils.”); 
return  1; 

} 

printf  ("What  is  the  file  name  for  color  map  output?  ") ; 
scanf  ("Xa" ,outf ilsname) ; 

if  ((file.ptl  ■  fopsn  (outfilename,  "sb"))  *■  TOLL) 

{ 

fprintf  (stdsrr,  "Cannot  open  output  file."); 
return  1; 

} 

printf  ("  What  is  the  file  name  for  data  file  output?  \n  " 

"  Use  a  sirtual  disk  address!!!  "); 
scanf  ("Xs",outfilename2) ; 

if  ((file  pt2  *  fopsn  (outf ilename2,  "sb"))  “  HULL) 

( 

fprintf  (stdsrr,  "Cazrnot  open  output  file."); 
return  1; 

> 

if  ((file  pt3  *  fopen  (outf ilenam*3,  "at"))  ■■  HULL) 

{ 

fprintf  (stdsrr,  "Cannot  open  input  file."); 
return  1; 

> 

//  Start  reading  the  header  information  in  the  first  three  lines  // 

//  See  if  first  sord  is  ’123456’  // 
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//  See  if  tha  data  is  compressed  // 

//  Sss  if  tha  color  naps  are  praaant  // 

//  Stora  tha  other  parts  of  the  sords  in  the  correct  parte  of  // 

//  tha  color  nap(s)  if  required  // 

//  // 


•hile  (!feof(file_pt)) 

{ 

//  Read  first  three  lines  // 

freed  (header .words, 1. 48, f ile.pt ) ; 

clrscrO ; 

array_indl-0; 

for  (line.ct  ■  0;  line.ct  <  3;  line_ct++) 

{ 

for  (loop_ct*0;  loop_ct<16;  loop_ct++) 

{ 

array.indl  »  array_indl++; 

} 

> 

//  Function  Call  // 

//  Decode  first  word  to  be  sura  it  counts  1234S6  // 

//  i.e.  A72E  -  lOlOIOllllOOlOllllO  -  __i  1 010 1  Oil  1 1001 101 1  UO  // 

//  -  1,2, 3, 4, 5, 6  // 

head[0]  «  dacode.header (header. words, 0); 
fprintf  (stdout ,"haad[0]  ■  Xx  \n  ",head[0]  ); 

(head [0]  !-  0xi72£) 

{ 

fprintf  (stderr,  "File  does  not  have  correct  header  word."); 
return  1; 

} 

//  Decode  second  aord  // 

headCl]  *  dacode.header (header. words, 1) ; 
fprintf  (stdout ."headCl]  »  Xx  \n  ",head[l]  ); 

//  Decode  third  word,  beginning  x  coord  // 

head [2]  ”  dacode.header (header. words, 2) ; 
fprintf  (stdout, "headC2]  ■  Xx  \n  " ,headC2]  ); 

//  Decode  fourth  word,  beginning  j  coord  // 

head [3]  »  dacode.header (header. words, 3) ; 
fprintf  (stdout, "head [3]  *  Xx  \n  ".head [3]  ); 

//  Decode  fifth  word,  ending  x  coord  // 

head [4]  *  dacode.header (header .words, 4); 
fprintf  (stdout, "head [4]  ■  Xx  \n  ",head[4]  ); 

//  Decode  sixth  word,  ending  y  coord  // 

head [5]  *  decode .header (header .words, 5) ; 
fprintf  (stdout, "head[5]  -  Xx  \n  ",head[5]  ); 
conp.data  ■  head[l]  k  0x80; 
reg_data  "  head[l]  k  0x40; 

■ap.l  ”  head[l]  k  0x10; 

■ap_2  “  headCl]  k  0x20; 
gflag  “  headCl]  k  0x08; 

fprintf  (stdout, "Xx\t  XxVt  Xx\t  Xx\t  Xx  \n",gflag,nap_l,nap_2, 
reg.data,  conp.data); 

if  (gflag  “  0x08)  //  File  contains  graphics  infornation  // 

{ 

if  (nap.l  ”  0x10)  //  File  has  top  color  nap  // 

{ 

store.nap.lO ; 

> 

if  (nap_2  “  0x20)  //  File  has  botton  color  nap  // 

store_nap_2(); 

} 
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if  (reg.data  "  0x40)  //  Data  ia  atorad  in  linaar  ordar  // 

{ 

straight. writaO ; 

> 

if  (coap.data  "  0x80)  //  Data  ia  atorad  in  conpraaaad  format  // 

{ 

unpack () ; 

} 

} 

} 

fcloaa  (fila.pt); 
fcloaa  (f ila.pt 1); 
fcloaa  (fila.pt2); 

> 

alaa 

{ 

printf  ("  What  ia  tha  fila  naaa  for  data  fila  output?  \n  " 

"  Uaa  a  virtual  diak  addrass!!!  "); 
acanf  ("Xs",outf ilenane2) ; 

} 

if  (Cfila_pt2  ■  fopan  (outf ilanana2,  "rb"))  “  HULL) 

f printf  (atdarr,  "Cannot  opan  output  fila."); 
raturn  1; 

} 

printf  ("  What  is  tha  fila  nana  for  corralation  data  fila  output?  \n  " 

"  Uaa  a  virtual  diak  addraaa!!!  "); 
acanf  ( "la" , outf ilsnana3) ; 

if  «fila_pt3  -  fopan  (outf ilanana3.  "at"))  —  TOLL) 

f printf  (atdarr,  "Cannot  opan  output  fila.”); 
raturn  1; 

} 

printf  ("  What  is  tha  fila  naaa  for  block  data  rectangles?  \n  " 

"  Uaa  a  virtual  diak  addraaa ! ! !  ") ; 
acanf  ("Xs",filenaae3); 

if  ((fila.pt 4  ■  fopan  (filanaaa3,  "rt"))  “  BULL) 

{ 

f printf  (atdarr,  "Cannot  opan  input  fila.”); 
raturn  1; 

} 

f printf (atdout,"  Do  you  wish  to  saa  tha  plot  in  ona  of  tha 
following  foraata?\n"); 

f printf (atdout,"  Linaar  Gray  Scale  (A)\nM); 

f printf (atdout,"  Log  Gray  Scale  (B)\nM); 

f printf (atdout,"  Color  (C)"); 

gats(inbuf) ; 
gata(inbuf); 

sscanf  (inbuf , "lie", tans); 

//  Exit  if  not  EGA  or  VGA  // 

//  Find  out  if  they  have  what  it  takas  // 
if  (set.gr aph()  »■  1) 

{ 

printf  ("This  progran  requires  EGA  or  VGA  graphical") ; 
axit(0) ; 

} 

cal c .coords O ;  //  Scale  to  graphics  resolution  in  use  // 

draw.acanO;  //  Draw  the  scan  fron  tha  fila  // 

draw.ractQ ;  //  Draw  tha  rectangles  around  data  // 
f printf (atdout,"  Find; (a)  \n"); 
f printf (atdout,"  average  \n"); 


f print 1 (stdout , "  (b)cor  1  \n"); 
gets (inbuf) ; 

■•cant  (inbuf ,"Xlc" .tansl) ; 

Result  ■  gsousa .Mreeet  O ; 

if (^(ansl  —  'a')  ||  (anal  —  *A’)II  (ansi  —  *b*)ll  (ansi  —  *B*)) 

if (  Rssult->pr«ssnt  } 

{ 

f print!  (stdprn,"  Oats  Fils  lass:  Xs  \n  ",outfilsnaas2); 

fprintf  (stdprn,"  Rectangle  File  Mane:  Xs  \n  " ,f ilenase3) ; 

f print!  (stdprn,"  Correlation  File  Rase:  Xs  \n  ".out filenames) ; 

log-gray (pal) ; 

aet  wit  anode  (  COPY.PUT  ); 

gnouae. Set .Cursor (  CIRCLE  ); 

gnouse.Nsbo*(  FALSE  ); 

Boil t an (  S78,  20,  SO,  20,  "Pt  1"  ); 

BozIten(  578,  50,  50,  20,  "Pt  2“  ); 

BozIten(  578,  80,  50,  20,  "Exit"  ); 
do 

( 

do  //  while  ansi  ■  a  or  A  // 

( 

if  ((ansi  —  ’a’)  II  (ansi— ’A')) 

{ 

surf ace. level () j  //  Pick  a  region  of  the  display  // 
statisticsO; 

fprintf  (stdprn."  AGAIR?  A  or  R  \n"); 

gataC  inbuf  ); 

sscanf  (inbuf ,"Xlc",Ransl) ; 

} 

>while(  (ansi  —  'a')  II  (ansi  —  ’A’)); 
do  //  while  ansi  -  b  or  B  // 

{ 

corral. lineO ;  //  Pick  a  region  of  the  display  // 

corral (); 

fprintf (stdprn,"  AGAIR?  B  or  I  \n"); 

gets(  inbuf  ); 

sscanf  (inbuf ,"Xlc",tansl); 

} while (  (ansi  —  ’b’)  II  (ansi  —  'B’)); 

}while(  (anal  —  ’a')  II  (ansi  —  ’A')  II 
(ansi  —  *b«)  ||  (ansi  —  *B»)); 

} 

} 

get.keyO;  //  Display  sassage  and  wait  for  key  press  // 

closegraphO ;  //  Close  graphics  system  // 

f close  (file.pt2); 
f close  (file.ptS); 
f close  (file.pt4); 
return  0; 

>  //  End  of  main  // 

//•••••***********•*••*•«**•*••••*•****•********••*•*•••**•****••***// 
unsigned  int  decode Jieader (unsigned  char  w[48],  unsigned  int  word.count) 

{ 

unsigned  char  header .word  ■  OzFO; 

unsigned  int  parti,  part 2,  part3,  part4,  il,  i2,  i3,  i4; 

//  The  indexes  are  reversed  because  the  segments  of  words  are  placed  in  the  U 
ft  header  in  reversed  order  // 

11  ■  word.count  *8+7; 

12  ■  word.count  *  8  ♦  5; 

13  ■  word.count  *8  +  3; 

14  ■  word.count  *8+1; 
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parti  ■  (w[il]  ft  header.vord)  «8 ; 
part  2  “  Cw[i2]  t  header.vord) <<4 ; 
part 3  ■  (w[i3]  ft  header.vord) ; 
part4  •  (w[i4]  ft  header.vord) »4 ; 
return  (parti  I  part 2  |  part 3  I  part4); 


//MM**»«*MM**M**mMM*«M**M*.*MMM**MMM*MM«*M*MM.// 


void  store.nap.l(void) 

{ 

unsigned  char  header .naskl  ■  OxOF; 
unsigned  char  header _aask2  ■  OxFO; 
unsigned  char  parti,  part2,  part 3 ; 
unsigned  long  part A,  partB,  partC; 
unsigned  int  il; 

//  Enter  the  data  fros  the  fir  t  48  lover  words  into  the  // 

//  first  24  color  sap  words  / 
for  (il  •  0;  il  <  24;  iH-fj 
< 

partB  “  (header _aask2  ft  header.vords [il*2]) ;  //  Lower  Qreen  // 
partC  “  (header _saskl  ft  header .words [i 1*2]) ;  //  Lower  Blue  // 
part A  ■  (header .saskl  ft  header. words [il*2+lj ) ;  //  Lower  Red  // 
color_sap_l[il]  »((partA«16)  |  (partB«4)  I  partC); 

} 

//  Get  the  rest  of  the  612  lower  words  and  convert  then  to  the  // 

//  loser  266  color  nap  1  values  // 
revind(f  ile_pt) ; 
fread  (aap_full,l,612,file_pt) ; 
array.indl“0; 

for  (line.ct  ■  0;  line.ct  <  32;  line.ct-*-*’) 

{ 

for  (loop.ct-O;  loop_ct<16;  loop_ct++) 

{ 

array _indl  ■  array_indl++ ; 

> 

} 

for  (il  -  24;  il  <  266;  il++> 

{ 

partB  a  (header _nask2  ft  nap.full[il*2]) ;  //  Lover  Green  // 
partC  ■  (header .naskl  ft  nap_full[il*2]) ;  //  Lover  Blue  // 
part A  ■  (header .naskl  ft  nap_full[il»2+ll) ;  //  Lover  Red  // 

color  .nap.  1  [il]  a((partA«l6)  I  (partB«4)  |  partC); 

} 

//  Get  the  512  upper  words  and  convert  then  to  the  // 

//  upper  256  color  nap  1  values  // 
fread  (nap_full,l,612,file_pt) ; 
array_indl-0; 

for  (line.ct  ■  0;  line.ct  <  32;  line.ctt-O 

{ 

for  (loop_ct“0;  loop_ct<16;  loop.ct-M') 

{ 

array.indl  ■  array_indl++; 

} 

} 

for  (il  -  0;  il  <  256;  il++) 

{ 

partB  ■  (header _nask2  ft  nap.full[il*2] ) ;  //  Upper  Green  // 

partC  •  (header .naskl  ft  nap.full[il«2]) ;  II  Upper  Blue  // 

part A  a  (header.naskl  ft  nap.full[il*2+l]) ;  //  Upper  Red  // 

color _nap_l[il]  -((partA«20)  I  (partB«8)  I  (partC«4)  I  color.nap.l [il] ) ; 

fprintf  (f ile.pt  1 , "Xc" , (Orff 000000  ft  color_nap_l[il])»24) ; 
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fpriatf  (file_ptl,"Xc", (0x00110000  ft  color_nap_l[il])»16) ; 
f print!  (file_ptl,"Xc",(OxOOOOffOO  ft  color_sap_l[il])»8) ; 
fprintl  (file_ptl."Xc",(OxOOOOOOff  ft  color.sap.lCil])) ; 


void  at ore.sap.2 (void) 

{ 

unsigned  char  header .saskl  ■  OxOF; 
unsigned  cJiar  hssdsr_nssk2  -  OxFO; 
unsigned  char  parti,  part2 ,  part 3 ; 
unsigned  long  parti,  partB ,  partC; 
unsigned  int  il; 

//  Oet  the  612  lover  words  and  convert  then  to  the  // 

//  lover  266  color  nap  2  values  // 
freed  (nap.full.l, 512, file.pt) ; 
array.indl*0; 

for  (line.ct  -  0;  line.ct  <  32;  line.ct*-*) 

{ 

for  (loop.ct-O;  loop_ct<16;  loop.ct*-*-) 

i 

array.indl  ■  array.indl*-*-; 

> 

> 

for  Cil  -  0;  il  <  266;  il-**-) 

{ 

partB  -  (header _sask2  ft  map_full[il*2] ) ;  //  Lover  Qreen  // 
partC  -  (header _naskl  ft  aap_fullCil*2]) ;  //  Lover  Blue  // 
parti  -  (header _naskl  ft  *ap_full[il*2+l]> ;  //  Lover  Red  // 

color _nap_2 [il]  *((parti«l6)  I  (partB«4)  I  partC); 

//  Oet  the  612  upper  words  and  convert  then  to  the  // 

//  upper  256  color  nap  2  values  // 
freed  (nap.full, 1,512,1 ile.pt) ; 
array.indl *0; 

for  (line.ct  ■  0;  line.ct  <  32;  line.ct*"*-) 

{ 

for  (loop_ct*0;  loop.ct<16;  loop.ct*-*-) 

{ 

array.indl  ■  array.indl*-*-; 

> 

> 

lor  (il  "0;  il  <  256;  il*-*-) 

{ 

partB  »  (header _nask2  ft  sap_full[il*2] ) ;  //  Upper  Oreen  // 

partC  "  (header .mask 1  ft  nap.full [il*2]  ) ;  //  Upper  Blue  // 

parti  »  (header .saskl  ft  nap.full  [il*2*-lj);  //  Upper  Red  // 

color jaap.2 Cil]  -((parti«20)  I  (partB«8)  |  (partC«4)  I  color.nap.2Cil]) ; 

f print!  (file_ptl,"Xc",  (0x11000000  ft  color_nap_2[il])»24); 

fprintl  (1 ile_ptl,"Xc",  (0x00110000  ft  color _nap_2[il])»16) ; 


f print!  (file.pt  1 , "Xc" ,  (0x00001100  ft  color.aap.2[il])»8) ; 
fprintl  (file_ptl,nXc", (0x00000011  ft  color_aap_2[il])) ; 

} 

> 

//  eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee// 
void  straight.vrite(void) 

{ 

//  Data  is  recorded  linearly  in  the  *. vch  file.  // 

//  Read  each  byte  and  write  the  byte  to  the  fila_pt2  file.  // 

//  The  data  vill  be  organized  in  512  lines  with  each  line  having  // 
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//  S12  data  bytea.  // 
array.  indl-O; 

for  (lina.ct  -  0;  lina.ct  <  512;  line  ct++) 

{ 

for  (loop_ct»0;  loop_ct<512;  loop_ct++) 

{ 

fraad  (ftdata.byte, aizeof (unsigned  char ),1, fila.pt); 
fprintf  (f ile_pt2 , "Xc" .data.byta) ; 

} 

} 

} 


//..•.....t..*,..*....*.,,,*..........,..*.*...*....*...............// 

void  unpack (Told) 

//  Data  la  coapraasad  in  tha  *.Ych  fila.  // 

//  Read  a  byte  and  decode  tha  inatruction  and  tha  count.  II 

//  If  tha  aord  aaya  nunber  continued  read  raat  of  nuaber  in  tha  next  byte.// 

//  If  tha  next  data  point  repreaanta  n  point*  write  n  point* .  // 

//  If  tha  next  data  point  repraaanta  tha  firat  of  n  point*,  // 

//  road  and  write  n  point a.  // 


//  Tha  data  will  be  organized  in  a  continuoua  lina.  // 

{ 

unaigned  long  total.count  *  0,  total.pointa; 

unaigned  char  aaq.taat  ■  0x80,  count .teat  "  0x40,  upper _aaak  •  Ox IF, 
lower _aaak  »  OxFF,  aaq.raault; 

unaigned  int  x.count  »  0,  y.count  -  0,  aaq.ct,  atraaa_ct; 
int  index; 

haad4  ■  head [4] + 1-head [2] ; 
heads  -  head[5]+l-haad[3] ; 
total.pointa  ■  head4  •  heads ; 

do 

i 

fraad  Chdata.byta , aizeof (unaigned  char) ,1, fila.pt); 
aaq.raault  *  data.byta  ft  aaq.taat; 

if  (aaq.raault  “  0x80)  //  A  lot  of  pointa  are  atorad  in  ona  •  // 

{ 

aaq.ct  •  upper _aaak  ft  data.byta;  //  t  of  pointa  // 
if  ((count.teat  ft  data_byte)  “  count .teat)  //  Count  cont’d?  II 

i 

II  Mora  count  in  next  word;  Read  it  and  ahift  tha  // 

//  higher  bita  to  tha  left  and  add  in  tha  aecond  // 

//  word  // 

fraad  (ftdata.byte,aizaof (unaigned  char) ,1, fila.pt) ; 
aaq.ct  ■  (aaq.ct  «  8)  I  (lowar.aaak  ft  data.byta); 

} 

fraad  (ftdata.byta, aizeof (unaigned  char) ,1, fila.pt); 

//  Store  that  data  word  in  aaq.ct  pointa  // 
for  (indax-0;  indax<aaq_ct ;  indax++) 

{ 

fprintf  (file_pt2,"Xc", data.byta); 

> 

total.count  +■  aaq.ct; 

} 

alaa  //  A  nuaber  of  different  waluaa  are  next  // 

atreaa_ct  ■  upper _aaak  ft  data.byta;  //  8  of  pointa  // 
if  ((count  teat  ft  data_byte)  ■■  count.teat)  //  Count  cont’d?  // 
{ 

//  Mora  count  in  nort  word;  Read  it  and  ahift  tha  // 

//  higher  bita  to  tha  left  and  add  in  tha  aecond  // 

//  word  // 
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fread  (tdata.byte.sizeof  (unsigned  char) ,  1,1  ile.pt) ; 
etreaa.ct  »  (streaa.ct  «  8)  t  (lover _aask  k  data.byte) ; 

} 

//  Stora  that  data  word  In  streaa.ct  points  H 
for  (index»0;  index<streaa_ct ;  index**) 

{ 

fraad  (kdata.byte.sizeof (unsigned  char) , 1 ,f il«_pt) ; 
f print!  (file_pt2,"Xcn .data.byte) ; 

} 

total.count  +■  atraaa.ct; 

} 

>  ahila  (total.count  <  total.pointa) ; 

/.count  ■  total.count  /  512; 
x.count  -  total.count  X  512; 

f print!  (stdout , "x.count  -  X3d  \t  /.count  ■  X3d" , x.count , /.count ) ; 

> 

int  set.graph(void) 

< 

int  graphdriver  -  DETECT,  graphsode,  error .code; 

//  Initialize  graphics  aystaa;  nust  be  EGA  or  VGA  // 
initgraph(tgraphdriver,  kgraphaode ,  "") ; 
error.code  ■  grapbresult () ; 
if  (error.code  !■  grOk) 

return(-l) ;  //  8o  graphics  hardware  found  // 

if  ((graphdriver  (•  EGA)  it  (graphdriver  I*  VGA)) 

{ 

closegraphO ; 
return  0; 

) 

retura(l);  //  Graphics  OK,  so  return  "true"  // 

} 

//aaaaaaaaaaaaaaaaaaeaaaaaeaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaeaaaaaaaaaaaaa// 

void  calc.coords (void) 

{ 

//  Set  global  variables  for  drawing  // 

sax.x  ■  getaaxxO;  //  Returns  saxiaua  x-coordinate  // 

aax.y  -  getaax/O ;  //  Returns  saxiaua  /-coordinate  // 

} 

//eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeaeeeeeeeeeeeeeaeeeeeeeeeeeeeeeeeeeeeee// 

void  draw.scan(void) 

{ 

//  initialize  graphics  and  local  variables  // 

//  of  the  setrgbpalette  function  // 

int  gdriver  ■  DETECT  ,  gaode ,  errorcode,  errorcodel; 
init  graph(kgdriver ,  A  gaode ,'"'); 
errorcode  ■  grapbresult (); 
if  (errorcode  !■  grOk) 

{ 

f print!  (stdout ,  "Graphics  error :  Xs\n" ,  grapherronsg  (errorcode)  ) ; 
fprintf (stdout ."Press  an/  key  to  halt:"); 
getchO ; 
ezit(l); 

} 


gatpalatta(tpal) ; 
getpalatte(*pal_old) ; 
pal.pt.old  -  gatdefaultpalattaO ; 

//  craata  gray  seal*  // 

if  ((ana  —  *a')  II  (ana  —  'A')) 
linaar .gray (pal) ; 
if  ((ana  —  *b')  II  (ana  —  ’B’)) 
l°g-gray (pal) ; 
nax.x  ■  gatnaxcolorO ; 
for^(acan_y  “  0:  acan_y  <  512;  acan.y**) 

fraad  (lina.acan.aixaof (unaignad  char) .512, fila_pt2) ; 
for  (»can_x»0;  acaa_z<512;  acan  x++) 

{ 

data.byta  »  lina_acan[acan_x] ; 
plot_paak  ■  paak.datact(data.byta) ; 
if (acan_y>31) 

{ 

putpixal  ( acan.x  ♦  64.  acan_y  -  32,  data_byta»4) ; 

} 

) 

//  Radiaplay  ???????  // 

f print! (atdout . "Paak-\nX3d\n" ,plot_paak) ; 

whila((ana— ’a')  1 1  (ana— ’bO  1 1  (ana— ’c’)  1 1  (ana— ’A’)  1 1  (ana—  *B*)  1 1 
(ana—’C’)  1 1  (ana— ’d’)  1 I  (ana— ’D')  I  i  (ana— ’  j ’)  1 1  (ana— ’ J’)  1 1 
(ana— ’r’)  1 1  (ana— ’R’)  1 1  (ana— ’z’)  1 1  (ana— ’Z’)) 

outtaxtxy(0,110.  "Redo?"); 
gata(inbnf); 

aacanf  (inbnf ,"Xlc",tana) ; 
if  ((ana  —  'a')  II  (ana  —  ’A’)) 
linear .gray (pal) ; 

alaa  if  ((ana  —  *b»)  II  (ana  —  *B’)) 
log-gray (pal) ; 

alaa  if  ((ana  —  ’c’)  II  (ana  —  *C’)) 

{ 

old_paletta(pal.old) ; 

) 

alaa  if  ((ana  —  *d’)  II  (ana  —  'D')) 

< 

naw.palatta(pal.old) ; 

} 

alaa  if  ((ana  —  *r ’)  II  (ana  —  ’R’)) 

draw.rectO ; 

> 

alaa  if  ((ana  —  ’j*)  II  (ana  —  *J’)) 

< 

joa.palatta(pal.old) ; 

} 

alaa  if  ((ana  —  •*')  II  (ana  —  'Z’)) 

< 

dra*_zooa(); 

} 

> 

} 

Toid  gat.kay(void) 

< 
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print*  (-  Prana  \n  " 

'•  any  kay  \n" 

"  to  or it  \nM) : 

gatchO ; 

} 

void  linaar_gr ay (struct  palattatypa  pal) 

{ 

int  i,  irsd ,  igraan,  iblua; 
satacti*apaga(0) ; 
for  (i-0;  i<pal .  sir  a;  i+-0 
< 

irad  “  i*4; 
igraan  ■  i«4; 
iblua  -  i*4; 

satrgbpalatta (pal. colors Ci] ,  irad.  igraan.  iblua); 

> 

} 

void  log_gr  ay  (struct  palattatypa  pal) 

< 

int  i.  irad.  igraan.  iblua; 
float  ifloat,  db.ratio; 

//  Calculata  a  log  acala  for  tha  gray  acala  // 
satacti*apaga(0) ; 

db.ratio  *  63.0  /  Iogl0(pal.oi*a+1) ; 
for  (i-0;  i<pal.siza;  i++) 

ifloat  -  i+1; 

irad  “  (db.ratio  a  loglO(ifloat)  ♦  0.6); 
igraan  *  irad; 
iblua  -  irad; 

satrgbpalatta(pal.colora[i] ,  irad,  igraan,  iblua); 

} 

} 

/ /aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa/ / 
int  paak.datact  (unaignad  char  data.point)  //  Find  bighaat  laval  in  plot  // 
{ 

if  (data_point>paak_to_now) 
paak_to.no*  -  data.point; 
raturn  paak_to.no* ; 


void  old_palatta (struct  palattatypa  pal) 

< 

satactirapaga(l) ; 
aatallpalatta(ftpal) ; 
satallpalatta(tpal) ; 

//  color  palatta  rad  groan  blua  // 

satrgbpalatta (pal. colors [0] ,  0  ,0  ,0); 
satrgbpalatta (pal. colors [l] ,  0,0,  63); 
satrgbpalatta  (pal.  colon  [3] ,  0  ,  63  ,  0); 
satrgbpalatta  (pal.  colon  [3] ,  0  ,  63  ,  63); 
satrgbpalatta  (pal .  colon  [4] ,  63  ,  0  ,  0); 
satrgbpalatta  (pal.  colon  [5] ,  63  ,  0  ,  63); 
satrgbpalatta(pal.colon[6] ,  44  ,  44  ,  44); 
satrgbpalatta(pal.colan[7] ,  44  ,  21  ,0); 
■atrgbpalatt*(pal.colon[8] ,  21  ,  21  ,  21); 
satrgbpalatta (pal. colors [9] ,  0,0,  32); 
satrgbpalatta  (pal.  colon  [10] ,  0  ,  32  ,  0); 
satrgbpalatta (pal. colors [1 1] ,  0  ,  32  ,  32); 
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•atrgbpalatta(pal.colors[12] ,  32  ,  0  .  0); 

•atrgbpalatta (pal. color* [13] ,  32  .  0  .  32) ; 

•atrgbpalatta (pal. colors [14]  ,  32  .  32  ,  0); 
•*trgbpal*tt*(pal.colors[lS] ,  63  ,  63  ,  63); 

} 

void  naa.palatta (struct  palettatypa  pal) 

{ 

•atactivapaga(l) ; 
satallpalatta(tpal) ; 

//  color  palatta  rad  groan  blu*  // 

•*trgbpal*tt«(pal. color* [0] ,  0  ,0  ,0); 

•atrgbpal*tt*(pal.colors[l] ,  32  ,  32  .  32); 

•*trgbpal*tt*(pal.colors[2] ,  63  ,  52  ,  63); 

•«trgbpalatta(pal.colors[3] ,  63  ,  0  ,  52); 
satrgbpalatta(pal.color*[4] ,  63  ,  0  ,  20); 
s*trgbpal*tt*(pal.color*[8],  63  ,  12  ,  0); 

•*trgbpal«tt*(pal.color*[6] ,  63  ,  44  ,  0); 

•atrgbpalatta (pal. colors [7] ,  44  ,  63  .0); 
satrgbpalatta (pal. colors [8] ,  12  ,  63  .0); 

•atrgbpalatta (pal. colors [9] ,  0  ,  63  ,  24); 

•atrgbpalatta(pal.colors[10] .  0  ,  63  ,  56); 
•atrgbpalatta(pal.colors[ll] ,  0  ,  32  ,  63); 
•atrgbpalatta(pal.colors[l2] ,  0,0,  63); 
satrgbpalatta(pal.colors[13] ,  32  ,  0  ,  63); 
•«trgbpalatta(pal.colars[14] ,  52  ,  0  ,  52); 
satrgbpalatt«(pal.colars[15] ,  20  ,  0  ,  20); 

} 

//**•****•**»********••**»******»»*****»***»**«**»**••********•**••**••// 

void  joa_palatta(stxuct  palattatypa  pal) 

{ 

satact ivapaga (1) ; 
satallpalatta(ftpal) ; 
satpalatta(0,21) ; 
aatpalatta(l,22); 
satpalatta(2,23) ; 

■atpalatta(3,24) ; 
satpalatta(4,25) ; 
satpalatta (5,26) ; 

■atpalatta(6,27) ; 
aatpalatta (7 , 28) ; 
satpalatta(8,29) ; 

■atpalatta(9,30); 
satpalatta(10,31) ; 
satpalatta (11 ,32) ; 
satpalatta(12,33) ; 
satpalatta ( 13 , 34) ; 
sat palatta (14, 35) ; 

■atpalatta(15,36) ; 

) 

//,,,..MMm****M*M*M*M*****MM.MM.MM*MM****MMM*.*// 

void  drav.zooa(Toid) 

[ 

int  dummy _x  ■  0,  dummy _y  “  0; 
outtartxy (0,140,  "X  START?”); 
gats(inbuf) ; 

•scant  (inbul , ”Xd” , Astart.z) ; 
outtartxy (0,170,  ”T  START?”); 
gats(iabui) ; 

•scan!  (inbul, ”Xd”, Astart.y) ; 


if  (start.!  <  0) 
start _z  ■  0; 
if  (start.y  <  0) 
start _y  “  0; 
if  (start.x  >”  256) 
start.!  ■  255; 
sod.!  -  start.!  ♦  256; 
if  (start.y  >■  256) 
start _y  ■  255; 
and.y  •  start.y  ♦  256; 
rssind(f  ila.pt  2) ; 

for  (scan_y  *  0;  scas.y  <  612;  scan_y«-0 

{ 

frsad  (lina.scan.sizaof (unsignad  char) ,612, fila_pt2) ; 
for  (scan.z»0;  acan_z<512;  scan_x++) 

{ 

data.byta  ■  lina_scan£scan_z] ; 
plot.psak  "  psak.dstsct (data.byta); 
if ( (scan_y>»atart_y) AA(*can_y<»and_y) ) 

{ 

if  (  (scan_z>«start_x)tt(acan_z<"and_x)  ) 

dussy.!  ■  2  •  (scan.!  -  start.!) ; 

dussy _y  »  2  *  (scan.y  -  start.y) ; 

putpizal(  (dussy _z  ♦  64)  ,  (dussy _y-32)  ,data_byts»4)  ; 

putpizaK  (dussy.!  ♦  65) ,  (dussy _y-32)  ,data_byts»4) ; 

putpizaK  (dussy.!  ♦  64) ,  (duasy_y-31)  ,data_byta»4)  ; 

putpizaK  (dussy.z  ♦  65)  ,  (dussy_y-31)  ,data_byta»4)  ; 

} 

> 

> 

} 

> 

//  *•»**•*******••*•**•*•*•***••*•****••••••••*•***•*••**•*•**•*•••••  // 

//  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa  // 
void  surf acs.lsvsl (void) 

{ 

//  Establish  ths  cursor  // 

Rasult  ■  gsousa.HrasatO ; 
satvritasoda(  C0PT.PUT  ); 
gsousa. Sat .Cursor (  CIRCLE  ); 
gsousa .Nshov(  FALSE  ); 
gsousa. Hshov(  TRDE  ); 
gsousa .Mxlinit (  64,  628  ); 

EXIT  -  FALSE; 

do  //  vhila(  EXIT  !-  TRUE  ) 

{ 

PositionR  ■  gsousa. Hprassad(  ButtonR  ); 

//  Cyds  through  salaction  until  satisfiad  // 
if(  TPos(  PositionR. zazis,  578,  628  )  ) 

{ 

if(  TPos(  PositionR. yazis,  20,  40  )  ) 
gat.pt 1(); 

if(  TPos(  PositionR. yazis,  50,  70  )  ) 
gat.pt  20; 

if(  TPos(  PositionR. yazis,  80,  100  )  ) 

EXIT  ■  TRUE; 

Draw  a  boz  around  tha  araa  // 
gsousa. Hsho«(  FALSE  ); 
if (uppar.laf t.z  >  lovar .right.!  ) 


// 
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{ 

teap.eoor  -  upper.left.x; 
upper.lef t.x  ■  lower.right.x; 
loser .right _x  -  teap.eoor; 

} 

if (upper.lef t.y  >  loser .right.y  ) 

{ 

teap.eoor  »  upper.lef t.y ; 
upper.lef t.y  ■  loser .right.y; 
lower.right.y  ■  teap.eoor; 

} 

rectangle (  upper.lef t.x,  upper.lef t.y , 
lower.right.x,  lower.right.y  ); 

} 

> 

}  while (  EXIT  !-  TRUE  ); 

} 

void  get_ptl(  void) 

//  Select  the  first  point  // 
gaouse .Hxl iait (  64,  575  ); 

do  //  while (PositionL. button_status**0  ) 

{ 

PositionL  -  gaouse . Hpreaeed (  ButtonL  ); 
if  (  PositionL. button_statua— 1  ) 

{ 

upper _l#ft_x  ■  PositionL. xsxis; 
upper.lef t.y  ■  PositionL. yex is ; 
gaouse . Kahow (  FiLSE  ); 
putpixel (upper. left_x,  upper _left_y,  0); 
gaous e. Kahow (  TRUE  ); 

} 

}shile(PositionL.button_status— 0  ); 
gaouse . Hzl iait (  64  ,  628  ); 

} 

//  eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 
void  get_pt2(  void) 

//  Select  the  second  point  // 
gaouse .Kxl iait (  64  ,  575  ); 

do  //  while  (Posit  ionL .  button.statua**  0) 

{ 

PositionL  ■  gaouse. Kpr eased (  ButtonL  ); 
if  (  PositionL. button_status**l  ) 

{ 

lower .right.x  -  PositionL. xsxis; 
lower.right.y  -  PositionL. y axis; 
gaouse . Hshow (  FALSE  ); 

putpixeKlower  .right  _x ,  lower  .right.y ,  0); 
gaouse. Hshow (  TRUE  ); 

} 

}  while  (Posit  ionL  .butt  on.st  at  us*»  0); 
gaouse. Mxl iait (  64  ,  628  ); 

} 

//••eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee. 

void  drsw.rect (woid) 

{ 

//  initialize  graphics  and  local  variables  // 

//  of  the  setrgbpalette  function  // 

freed  (An .boxes , sizeof (unsigned  char) , 3 , f ile_pt4) ; 
boxes*  atol(n.boxea) ; 


for  (n.box  “  0;  n.box  <  bozos;  n.box-*"*') 

{ 

froad  (ftxl.sizeof (unsigned  char) ,4,filo_pt4) ; 
froad  (kyl.sizeof (unsigned  char) ,4,f ile_pt4) ; 
tread  (ftx2,sizeof (unsigned  char) ,4,f ile_pt4) ; 
fread  (*y 2, sizeof (unsigned  char) ,4,f ile_pt4) ; 
upper.left.x  “  atoi(zl) ; 
upper.left.y  ■  atoi(yl); 
lower .rigbt.x  ■  atoi(z2); 
loser .right.;  ■  atoi(y2); 
rectangle (  upper_left.x,  upper.left.y, 
loser _right_x,  loser_right_y  ); 

} 

> 

soid  statistics(roid) 

//  Rewind  the  data  file  and  find  the  awerage  value  of  the  level  // 

{ 

fpriatf  (stdprn, "  XI  -  Xd  Y1  -  Xd  X2  -  Xd  T2  «  Xd  \n  " , upper.lef t_ 

upper.left.y,  loserjright.z,  loser .right.;  ); 
zaaz.liait  »  loser .right _z  -  64; 
ysaz.liait  -  lower .right.;  ♦  32; 
average  “  0; 
average_2  *  0; 
st .dev  “  0; 
n.points  »  0; 
n.points.off  -  0; 
plot .peak  -  0; 
peak.to.nos  «  0; 
rewind  (file_pt2) ; 

for  (scan.;  -  0;  scan.;  <■  ymaz. limit;  sean.ye-O 

i 

fread  (lino.scan , sizeof (unsigned  char) , 612, f ile_pt2) ; 

if(  (scan_y>“ (upper.lef  t_y*-32) )  ft*  (scan_y<» (lower .right _y+32))  ) 

{ 

for  (scan_z-0;  scan_z<“zmaz_li»it ;  scan_x++) 

if (  ( sc an_x>“ (upper .left.x  -  64))  k* 

(scan.xO (loser .right _z  -  64))  ) 

{ 

data  .byte  ■  line  .scan  [scan,  x]  ; 

data.float  ■  data. byte; 

plot.peak  ”  peak.detect(data.byte) ; 

average  *  average  ♦  data.float; 

average_2  *  average.2  ♦  (data.float  *  data.float); 

n.points  "  n.points  ♦  1; 

> 

) 

} 

} 

average  »  average  /  n.points; 

average. 2  -  (  average_2  -  n.points  e  average  e  average  ) ; 

//  Find  the  standard  deviation  // 

if (  average_2  <«  0  ) 

{ 

st.dev  -  0.0; 

} 

else 

{ 

st.dev  ■  sqrt  (  (  average_2  /  n.points  )  ) ; 

} 


//  Print  thin  intonation  out  // 

fprintf  (atdprn,"  I unbar  of  Point*  ■  %6.0f  \n  " ,n_point»  ); 

fprintl  <»tdprn,"  Peak  ■  Xd  \n  ".plot.peak  >; 

fprintf  (atdprn , "  Moan  -  X6.21  \n  " .average  ); 

fprintl  (atdprn,"  Standard  Deviation  “  X6.21  \n  *, at.dev  ); 

} 

void  corral (void) 

{ 

//  Rewind  the  data  file  and  find  the  correlation  length  of  the  data  // 
fprintf  (atdprn, "  XI  ■  Xd  Yi  ■  Xd  X2  *  Xd  Y2  ■  Xd  \n  ",upper_l*ft 

upper.laft.y,  lover _right_x ,  lover .right.y  ); 
xnax.linit  -  lover .right.*  -  64; 
ynaz.linit  ■  lover .right.y  ♦  32; 
for  (i.a*  0;  i_a<*  S12;  i_a++) 

{ 

dat  a_ array  [i_a]  “0.0; 

) 

average  •  0; 
average. 2  “  0; 
at.dev  ■  0; 
npoint  ■  0; 
n.pointa  -  0; 
n.pointa.off  ■  0; 
plot .peak  ■  0; 
peak_to.no*  ■  0; 
r«vind(f ile_pt2) ; 

for  (acan.y  -  0;  acan.y  <-  ynaz.linit;  acan.yM-) 

{ 

freed  (lina.acan.aizeof (unaigned  char) .512, file.pt 2) ; 

if (  (acan_y>-(upper_left_y+32))  U  (*can_y<-(lo*er_right_y+32))  ) 

{ 

for  (scan.z-0;  *can_z<“xnaz_linit ;  acan_z++) 

{ 

if (  (»can_x>» (upper .left.x  -  64))  U 
(acan_z<» (lover .right _z  -  64))  ) 

{ 

data. byte  »  line_acan[*can_x]  ; 
dat a. float  -  data.byte; 
plot .peak  ■  peak.det wet (data.byte) ; 
average  *  average  ♦  data.float; 

average_2  ■  average. 2  *  (data.float  •  data.float); 
n.pointa  ■  n.pointa  ♦  1; 
data.array [npoint]  »  data.float; 
npoint  -  npoint  +  1; 

> 

} 

} 

> 

average  »  average  /  n.pointa; 

average_2  “  (  average.2  -  n.pointa  *  average  *  average  ); 

//  Find  the  standard  deviation  // 
if(  average.2  <■  0  ) 

at.dev  •  0.0; 

} 

else 

{ 

at.dev  ■  sqrt(  (  average.2  /  n.pointa  )  ); 

} 

//  Find  the  correlation  length  // 

//  First  repeat  the  array  in  data.array 
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lor  (i.a"  0;  i_a<“  npoint-1;  i_a++) 

{ 

data.array  [i_a]  ■  data.array  [i.a]  -  ayarage ; 
data. array [i_a  ♦  npoint]  ■  data.array [i.a] j 
} 

lor  (i_a“  0;  i_a<“  npoint;  i_a++) 

{ 

lor  (  j_a«  0;j_a  <■  npoint;  j_a++) 

{ 

corr array [i.a]  »  corrarray[i_a]  ♦ 
data.array [j.a]  •  data.array [j.a  +  i.a]; 

} 

corr array [i.a]  ■  corrarray [i.a]  /  n.points; 

> 

lor  (j.a*  0;j_a  <-  npoint;  j.a-t-O 

{ 

f.data  ■  corrarray  [j.a]  ; 

Iprinti (f ile_pt3 ,  "XI  \n",  f.data) ; 

} 

cor.yalue  ■  corrarray [0]  *  0.367879;  //  a*-l 

i.a  ■  0; 

do  //  while  corralation  Talua  im  greater  than  cor .value  // 

{ 

i.a  «  i.a  +  1; 

}*hila(  (corrarray [i.a]  >  cor_valua)kl(i_a  <  npoint)); 
indax.cor  «  i.a; 

//  Print  thia  inlornation  out  // 

fprintf  (atdprn,"  limber  of  Points  *  X6.0f  \n  ”, n.points  ); 

fprintf  (atdprn,"  Peak  *  Xd  \n  ", plot .peak  ); 

fprintf  (stdprn,"  Mean  “  X6.2f  \n  ".average  ); 

fprintf  (stdprn."  Standard  Deviation  *  X6.2f  \n  ",st_dev  ); 

fprintf  (stdprn,"  First  Corralation  Valua  "  X6.2f  \n  ", corrarray [0]  ); 

fprintf  (stdprn,"  Corralation  Langth  Valua  *  X6.2f  \n  ", 

corrarray [indax.cor]  ) ; 

fprintf  (stdprn,"  Corralation  Length  in  Pixels  “  Xd\n  ".indax.cor); 

> 

//  eeeeeeeeeeeeeeeeeeeeeeeeeeeeeaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa  // 

roid  correl.line(Toid) 

{ 

//  Establish  the  cursor  // 

Result  *  gnouse.  HresetO; 
seter it anode (  C0PT.PUT  ); 
gnouse. Set. Cursor (  CIRCLE  ); 
gnouse . Mshow (  FALSE  ); 
gnouse. Hshov(  TRUE  ); 
gnouse. Mxlinit(  64,  628  ); 

EXIT1  *  FALSE; 

do  //  while (  EXIT1  !-  TRUE  ) 

{ 

PositionR  ■  gnouse. Mpressed(  ButtonR  ); 

//  Cycle  through  selection  until  satisfied  // 
if(  TPos(  PositionR.xaxis,  578,  628  )  ) 

if(  TPos(  PositionR. yaxis,  20,  40  )  ) 
get.ptlO; 

if(  TPos(  PositionR. yaxis,  50,  70  )  ) 
get_pt2() ; 

if(  TPos(  PositionR. yaxis,  80,  100  )  ) 

{ 

EXIT1  -  TRUE; 

Draw  a  box  around  the  area  // 


// 
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g»ou»«.I!*how(  FALSE  ); 

//  Find  til*  long  dinansion  // 

d#1-x  *  losar.right.x  -  uppar.laft.x  ♦  i; 
d#1-y  “  lo»«r.right_y  -  upp*r_l«ft  y  1; 
if  (dal_y  >-  dal_x) 

{ 

l«»«r.right_x  ■  uppar.laft.x; 
dal.x  -  lj 
> 

■Isa 

{ 

lo»«r_right_y  "  uppax.laft.y; 
dal.y  »  1; 

} 

ractanglaC  uppar.laft.x,  upp*r_l*ft_y, 
loaar.right.x,  lowar.right.y  ); 

} 

}  whila(  EIIT1  !-  TRUE  ); 


